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The  first  step  in  our  program  was  to  demonstrate  that  the  magnetically  compen¬ 
sated  supersonic  beam  was  a  viable  concept,  by  achieving  a  very  large  single  pass 
absorption  coefficient  and  narrow  linewidth  in  a  weak  atomic  transition.  Briefly,  the 
compensated  beam  system  consists  of  a  tantalum  tube  oven,  filled  partially  with  the 
atomic  metal  of  interest,  and  heated  by  a  current  of  up  to  500  A.  Drilled  in  the  side  of 
the  tube  is  a  pinhole  of  50-150  microns  diameter  through  which  a  supersonic  atomic 
beam  emerges  after  adiabatic  expansion.  This  beam  has  the  narrow  spread  of  thermal 
speed  due  to  the  cooling  effect  of  the  expansion.  A  laser  beam  intersects  the  atomic 
beam  at  right  angles  to  the  central  axis  of  the  emerging  cone  of  atoms.  In  the  absence 
of  any  magnetic  compensation,  atoms  which  travel  at  different  angles  with  respect 
to  the  laser  beam  experience  different  Doppler  shifts.  In  this  case,  the  absorption 
linewidth  for  an  interaction  path  length  of  4.4  cm  is  about  500  MHz  for  a  typical  vis¬ 
ible  transition.  The  interaction  is  limited  to  those  atoms  which  are  Doppler  shifted 
within  a  natural  linewidth  (or  transit  width)  of  resonance.  In  order  to  dramatically 
increase  the  number  of  interacting  atoms,  a  pair  of  opposing  magnet  coils  is  placed 
on  the  laser  beam  axis  on  either  side  of  interaction  region.  For  the  atomic  transitions 
employed  in  the  experiments,  which  are  all  7  =  0  ground  state  to  7  =  l,iVf  =  1 
excited  state  optical  transitions,  the  magnet  coils  have  the  effect  of  producing  a  field 
gradient  which  shifts  the  excited  state  as  a  function  of  the  atomic  position  along  the 
laser  axis.  Since  the  Doppler  shift  varies  approximately  linearly  in  position  along 
the  l2iser  beam  axis,  the  strength,  of  the  magnetic  field  gradient  generated  by  the 
opposing  magnet  coils  can  be  adjusted  to  cancel  the  Doppler  shift  over  an  extended 
region  along  the  Iciser  beam  axis.  This  leads  to  a  very  narrow  resonance  line  and 
tremendously  enhanced  absorption  which  is  ideal  for  nonlinear  optics  experiments. 

This  method  wais  demonstrated  initially  using  the  ^5o  — Pi  transtion  in 
at  556  nm.  With  a  radiative  lifetime  of  875  ns,  the  absorption  coefficient  for  a  wezdc 
laser  beam  traversing  a  4.4  cm  path  of  atoms  was  so  large  that  it  did  not  go  below 
50%  until  the  laser  was  detuned  more  than  50  MHz  from  resonance. 

In  order  to  carefully  study  the  lineshape  for  the  magnetically  compensated  beam, 
it  was  realized  that  a  much  weaker  transition  was  needed.  For  this  purpose,  the 
*5o  — Pi  transition  in  ®®5r  at  690  nm  W2is  selected.  With  a  radiative  lifetime  of 
22  fisec,  this  transition  seemed  ideal  for  the  study,  and  opened  up  the  possibility  of 
coherent  interaction  of  the  atomic  sytem  with  the  la«er  field  over  extended  distance 
along  the  atomic  beam  axis.  Unfortunately,  operation  of  the  ring  dye  laser  at  this 
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wavelength  proved  to  be  a  formidable  task.  Pyridine  2  optics  were  used  with  Pyridine 
1  dye.  When  these  were  installed  in  the  699-21  dye  laser,  it  was  not  possible  to  obtain 
lasing.  This  resulted  in  a  substantial  time  loss.  Finally,  we  inserted  a  100%  reflector 
in  the  laser  output  port  and  converted  to  standing  wave  operation.  Once  the  initial 
alignment  Wcis  completed,  a  4%  output  coupler  was  inserted  and  500  mW  at  690  nm 
was  obtained  open  cavity  using  6  W  of  all  lines  argon  ion  pump  power.  By  using 
a  2%  output  coupler,  we  finally  obtained  100  mW  of  stable  single  mode  power  with 
only  5  W  of  pump  power,  which  enabled  us  to  begin  the  lineshape  experiments  in 
strontium. 

The  data  obtained  in  the  strontium  experiments  demonstrated  beautifully  the 
magnetically  compensated  beam  method  and  yielded  good  lineshapes  which  could  be 
compared  to  predictions.  An  absorption  enhancement  of  29:1  wm  obtained  with  the 
compensating  coils  to  achieve  a  50%  absorption  and  a  compensated  linewidth  of  18 
MHz.  The  atomic  density  in  the  4.4  cm  interaction  path  was  estimated  at  1.2  x  10". 
The  lineshapes  obtained  were  in  good  agreement  with  theoretical  predictions  which 
included  the  frequency  chirping  due  to  the  atomic  motion  across  the  magnetic  field 
gradient. 


Publications  based  on  this  Research 

1)  K.D.  Stokes,  C.  Schnurr,  J.  Gardner,  M.  Marable,  S.  Shaw,  M.  Goforth,  D.E. 
Holmgren,  and  J.  E.  Thomas,  “Magnetically  Compensated  Supersonic  Beams  for 
Nonlinear  Optics,”  Optics  Letters  14,  1326  (1989). 

II.  CONTINUOUS  SPATIAL  PHOTON  ECHOES 

In  demonstrating  the  magnetic  compensation  method  in  atomic  strontium,  it  was 
appreciated  that  the  long  radiative  lifetime  of  the  excited  electronic  state  opened 
up  the  possibility  of  achieving  continuous  coherent  interaction  of  the  atomic  system 
with  laser  fields  over  large  distances  along  the  atomic  beam  axis.  One  important 
application  is  the  development  of  continuous  spatial  photon  echoes  in  atomic  beams. 
Continuous  echoes  can  be  used  for  nonlinear  optical  processing,  including  spatial 
phase  conjugation  and  correlation  or  convolution  of  pairs  of  optical  fields.  In  addition, 
continuous  generation  of  dense  macroscopic  optical  coherence  in  atomic  beams  has 
applications  in  a  variety  of  fundamental  experiments,  surh  the  evolution  of  optical 
coherence  in  cavities,  and  the  study  of  optical  quantum  noise. 
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Due  to  the  difficulty  we  experienced  in  the  strontium  experiments,  we  sought  a 
different  atomic  system  for  the  echo  experiments.  The  appropriate  system  needed  a 
simple  two  level  atomic  transition  in  the  visible  region  with  a  somewhat  longer  spon¬ 
taneous  lifetime  than  ytterbium.  After  some  searching,  the  samarium  atomic  system 
presented  itself  as  a  good  choice.  For  the  first  experiments,  the  ^Fo  Di  transition 
at  599  nm  was  chosen.  With  a  1.7  fisec  radiative  lifetime,  it  permitted  coherence 
to  propagate  over  millimeter  distances  along  the  atomic  beam  axis.  Further,  the 
transition  permitted  the  use  of  R6G  dye,  greatly  simplifying  the  dye  laser  operation. 

In  the  experiments,  the  atomic  beam  is  magnetically  compensated  by  taking  ad¬ 
vantage  of  the  excited  state  Zeeman  tuning,  for  which  the  g-factor  is  2.83.  For  a  1 
cm  interaction  path  length  with  the  laser  fields,  the  Doppler  width  without  compen¬ 
sation,  100  MHz,  is  reduced  to  5  MHz.  The  atoms  cross  two  cw  focussed  laiser  beaun 
regions,  each  of  50  micron  diameter.  This  yields  a  transit  time  bandwidth  somewhat 
larger  than  the  compensated  Doppler  width.  In  this  way,  nearly  all  atoms  interact 
strongly  with  the  laser  fields.  The  field  intensities  were  adjusted  to  yield  nominally 
ir/2  and  ir  pulse  areas  in  the  atom  frame,  to  optimize  the  echo  amplitude.  With  the 
laser  fields  separated  by  0.5  mm,  a  strong  continuous  echo  field  was  observed  1  mm 
downstream  from  the  first  laser  field.  Unlike  previous  work  on  continuous  coherent 
radiation  by  atoms  which  traverse  separated  fields,  the  echo  fields  generated  in  the 
present  experiments  could  be  viewed  with  the  unaided  eye.  Prior  work  utilized  het¬ 
erodyne  beats  with  a  strong  local  oscillator  to  make  mezisurements.  A  permanent 
record  of  the  spatial  echo  signals  in  our  experiments  was  obtained  by  imaging  the 
echo  field  onto  a  Reticon  diode  array. 

The  measurements  yielded  a  somewhat  unexpected  nonexponentiaJ  decay  charac¬ 
teristic  for  the  spatial  echo  signal  as  a  function  of  separation  between  the  input  laser 
fields.  This  was  found  to  be  due  to  the  magnetic  compensation,  which  causes  the 
atomic  radiation  frequency  to  chirp  as  the  atoms  move  through  the  field  gradient.  A 
simple  theory  taking  the  chirp  effect  into  account  yielded  predictions  in  qualitative 
agreement  with  the  data. 

The  experiments  demonstrated  that  very  dense  optical  coherence  can  be  gener¬ 
ated  and  transported  downstream  in  a  magneitcally  compensated  beam.  This  method 
opens  up  the  possibility  of  a  variety  of  fundamental  and  applied  physics  experiments 
as  discussed  above. 


Publications  based  on  this  Research 

1)  C.  Schnurr.  K.D.  Stokes,  G.R.  Welch,  and  J.E.  Thomas,  “Continuous  Spatial  Pho- 
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ton  Echoes.”  Optics  Letters  15,  1098  (1990). 

2)  J.E.  Thomas,  K.D.  Stokes,  C.  Schnurr,  and  G.R.  Welch,  “  Continuous  Spatial 
Photon  Echoes,”  IQEC  ’90  (Anaheim,  May  1990). 

III.  PRECISION  ATOMIC  POSITION  MEASUREMENT 
USING  OPTICAL  FIELDS 

The  success  in  applying  magnetic  field  gradients  to  compensating  the  Doppler 
broadening  of  diverging  supersonic  beams  suggested  that  it  might  be  fruitful  to  pursue 
other  applications  of  field  gradients.  Motivated  by  the  lack  of  position  measurement 
techniques  for  moving  atoms  in  atomic  beams,  we  have  developed  a  novel  atomic 
imaging  technique  based  on  resonance  imaging  in  intense  field  gradients.  The  purpose 
of  the  methods  which  are  being  developed  is  to  replace  the  mechanical  hot  wire  and  slit 
methods  generally  used  for  measuring  atomic  beam  intensity  distributions  with  high 
resolution  optical  techniques.  It  is  expected  that  the  new  techniques,  which  ultimately 
can  achieve  uncertainty  limited  position  resolution  in  the  nanonmeter  region,  will  find 
wide  applications  in  the  newly  emerging  field  of  atomic  de  Broglie  optics. 

The  basic  method  is  use  a  spatially  varying  potential  to  make  the  resonance  fre¬ 
quency  of  simple  atoms  position  dependent,  analagous  to  magnetic  resonance  (MR) 
imaging.  Unlike  conventional  MR  imaging,  the  method  employs  Raman-  induced 
resonance  imaging  using  optical  fields.  The  use  of  optical  fields  permits  the  measure¬ 
ments  to  be  made  in  very  small  volumes  with  very  large  spatially  varying  potentials. 
This  leads  to  quantum  mechanical  limits  on  the  measurement  accuracy. 

In  the  initial  experiments,  an  atomic  samarium  beam  is  collimated  to  1  mr  and 
passes  through  three  regions: 

i)  An  optical  pumping  beam,  resonant  on  the  ^Fi  Fq  transition  at  571  nm.  This 
beam  empties  the  M  =  1  sublevel  of  the  "Fi  electronic  ground  fine  structure  state, 
leaving  population  in  the  M  =  Q  sublevel.  This  level  is  used  as  the  initial  state  for  a 
Raman  transition  back  to  the  M  =  1  sublevel. 

ii)  A  Raman  region  containing  two  optical  fields,  both  of  which  are  100  MHz  off 
resonant  with  the  excited  state.  This  region  induces  Raman  transitions  from 
the  M  =  0  initial  state  to  the  A/  =  1  final  state,  which  is  initially  empty.  In  this 
region,  a  100  G  uniform  magnetic  field  plus  a  500  G/cm  field  gradient  are  applied. 
The  net  effect  of  these  two  fields  is  to  shift  the  A/  =  1  ground  sublevel  at  a  rate 
of  10*  Hz/cm.  The  optical  Raman  fields  are  generated  by  acoustooptic  modulation 
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from  a  common  laser  field  to  eliminate  laser  jitter  and  achieve  transit  time  limited 
frequency  resolution  on  the  Raman  transition.  With  an  atomic  beam  speed  of  8.6 
xlO^  cm/sec  and  a  Raman  region  diameter  of  2.3  mm,  transit  limited  resolution  of 
230  kHz  (FWHM)  is  predicted  and  obtained  in  the  experiments  when  the  gradient 
magnet  is  off. 

With  the  gradient  magnet  on,  only  atoms  in  the  spatial  region  tuned  within  230 
kHz  of  resonance  with  the  Raman  fields  make  a  Raman  transition  from  the  initial 
M  =  0  state  to  the  final  M  =  1  state. 

iii)  An  optical  field  resonant  with  the  M  —  I  — Fq  interrogates  the  itf  =  1  final 
state.  The  number  of  atoms  in  the  final  state  is  proportional  to  the  number  of  atoms 
in  the  initial  state  near  the  spatial  region  where  resonance  with  the  Raman  fields 
occurs.  By  measuring  resonance  fluorescence  intensity  in  this  region  eis  a  function  of 
the  uniform  magnitic  field  strength,  the  spatial  region  resonant  with  the  Raman  fields 
is  swept  across  the  atomic  distribution,  yielding  the  atomic  density  as  a  function  of 
position  across  the  atomic  beam. 

At  present,  spatial  resolution  has  been  studied  by  measuring  the  shadow  of  a  25 
micron  diameter  gold  coated  tungsten  wire  placed  just  before  the  Raman  fields.  To 
date,  the  best  resolution  is  about  2.5  microns,  which  is  close  to  the  theoretical  limit 
of  1.7  microns  due  to  the  velocity  of  the  moving  atoms  along  the  measument  axis. 

Much  higher  resolution  will  be  obtained  with  the  next  generation  systems  which 
employ  large  spatially  varying  light  shifts  instead  of  magnetic  field  gradients.  This 
method  permits  final  state  tuning  rates  of  10**  Hz/cm  over  micron  regions.  Since 
the  initial  state  is  not  shifted  by  the  applied  spatially  varying  potential,  momentum 
is  imparted  to  the  atom  only  after  it  makes  a  transition  into  the  final  state.  In  this 
C2ise,  for  adequately  collimated  atomic  beams  (50  fix  half  angle)  it  is  easy  to  show 
that  the  momentum  spread  imparted  by  the  force  due  to  the  field  gradient  and  the 
position  resolution  obey  the  Heisenberg  uncertainty  relation.  In  fact,  for  a  Raman 
region  comprised  of  laser  beams  with  a  gaussian  spatial  intensity  distribution,  it  is 
possible  to  put  atoms  into  the  final  state  in  a  gaussian  minimum  uncertainty  wave 
packet. 


Publications  based  on  this  Research 


1)  J.E.  Thomas,  “Uncertainty-Limited  Position  Measurement  of  Moving  Atoms  using 
Optical  Fields,”  Optics  Letters  14,  1186  (1989). 


2)  J.E.  Thomas,  “Quantum  Theory  of  Atomic  Position  Measurement  using  Optical 
Fields,”  Phys.  Rev.  A42,  5652  (1990). 

3)  J.E.  Thomas,  K.D.  Stokes,  C.  Schnurr,  J.  Gardner,  and  G.R.  Welch,  “Atomic  Po¬ 
sition  Measurement  and  Localization  using  Optical  Fields,”  De  Broglie  Optics  Work¬ 
shop,  (Santa  Fe,  N.M.  1990). 

4)  G.R.  Welch,  K.D.  Stokes,  C.  Schnurr,  J.  Gardner,  an  J.E  Thomas,  “Uncertainty 
Limited  Atomic  Position  Measurement  using  Optical  Fields,”  Symposium  on  Atomic 
Optics,  Annual  Meeting  of  the  Optical  Society  of  America,  Boston,  MA  (November, 
1990). 


Report  of  Inventions  based  on  this  Research 

1)  Precision  Optical  Measuring  System  for  Atoms  in  Beams 
(submitted  to  Duke  office  of  Technology  Transfer,  file  ^00635) 
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Magnetically  compensated  supersonic  beams 
for  nonlinear  optics 
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Doppler  frequently  shifts  for  a  laser  field  interacting  with  a  diverging  supersonic  atomic  beam  are  canceled  using 
spatially  varying  Zeeman  shifts.  Dense  atomic  beams  with  long  interaction  path  lengths  and  narrow  linewidths  are 
obtained  for  spectroscopic  and  nonlinear-optics  applications. 


It  is  well  known  that  supersonic  beams  are  useful  for 
producing  dense  sources  of  atoms  and  molecules  for 
spectroscopic  applications.^-^  In  such  beams,  adia¬ 
batic  expansion  causes  the  temperature  to  drop  so  that 
the  velocity  distribution  relative  to  the  average  stream 
speed  is  narrow.*-^  This  low-temperature  characteris¬ 
tic  has  been  exploited  in  a  variety  of  experiments  to 
produce  weakly  bound  species,  such  as  van  der  Waals’s 
molecules,^  and  to  produce  clusters.^  For  optical  laser 
spectroscopy,  however,  the  wide  divergence  angle  of 
the  expanding  beam  leads  to  substantial  Doppler 
broadening,  which  is  reduced  in  many  experiments  by 
collimation.  Collimation  is  adequate  whenever  the 
signal-to-background  ratio  or  experimental  require¬ 
ments  permit  a  short  interaction  path  length.  In 
many  experiments  this  is  not  the  case,  and  alternative 
means  of  reducing  the  linewidth  have  been  developed. 
Recently,  pulsed  slit  nozzles  have  been  used  to  reduce 
Doppler  broadening  of  a  long-path-length  (2.5-cm) 
molecular  beain  by  exploiting  the  fact  that  the  stream¬ 
lines  of  the  beam  tend  to  be  perpendicular  to  the  slit 
except  at  the  edges.^  Linewidth  reductions  of  three¬ 
fold  to  sixfold  in  infrared  spectra  have  been  obtained 
in  this  way.^  While  this  technique  is  general,  the  slit 
conflguration  requires  pulsed  operation  to  reduce 
source  consumption  and  pumping  requirements. 
Furthermore,  the  linewidth  reduction  is  modest  for 
some  applications. 

In  nonlinear-optics  applications,  it  is  often  impor¬ 
tant  to  have  a  high-density  beam  and  a  long  interac¬ 
tion  path  with  a  narrow  linewidth  in  the  optical  region. 
In  this  Letter  a  new  technique  is  demonstrated  that 
achieves  this  goal  by  using  a  spatially  varying  Zeeman 
shift  to  compensate  the  Doppler  broadening  of  a  di¬ 
verging  supersonic  beam.^  A  greater  than  20-fold 
linewidth  reduction  and  peak  absorption  enhance¬ 
ment  are  obtained  on  the  *So  ®Pi  transition  at  690 
nm  in  *®Sr,  for  which  the  radiative  lifetime  is  22.4  nsec. 
Absorptions  of  greater  than  50%  and  a  linewidth  of 
~18  MHz  are  obtained  with  a  4.4-cm  interaction  path 
length. 

The  *So  -*  *Pi  transition  at  690  nm  in  ®®Sr  was 
chosen  for  the  demonstration  because  of  its  long  22.4- 
/isec  radiative  lifetime  and  correspondingly  weak  tran¬ 
sition  strength.  In  this  way,  the  compensated  line 


shape  can  be  accurately  studied  at  high  density  and 
long  path  length  without  the  sample’s  becoming  opti¬ 
cally  thick. 

liie  basic  features  of  the  compensated  beam  system 
are  shown  in  Fig.  1.  A  tantalum  oven  is  heated  with  a 
few  hundred  amperes  of  ac  current  to  obtain  a  high 
density  of  active  species.  A  supersonic  beam  emerges 
from  a  ISO-^m  pinhole  nozzle  drilled  in  the  side  of  the 
tantalum  tube  and  passes  through  a  2.5-cm  aperture 
located  8  cm  downstream  from  the  nozzle.  The  corre¬ 
sponding  divergence  half-angle  is  0.156  rad.  The  spe¬ 
cies  chosen  for  study,  Sr  and  Yb,  have  two  valence 
electrons  and  negligible  dimer  formation.  Hence,  no 
carrier  gas  is  needed  for  the  experiments,  and,  since 
the  met^  are  condensable,  only  a  small,  5-cm  (2-in.) 
diffusion  pump  is  used  for  the  vacuum  S3rstem.  A 
circularly  polarized  laser  beam  of  ~0.5-mm  1/e  inten¬ 
sity  radius  crosses  the  atomic  beam  14  cm  from  the 
nozzle  so  that  the  interaction  path  length  is  4.4  cm.  In 
this  case,  taking  the  laser  beam  propagation  direction 
as  the  quantization  i  axis,  the  c/*0,M  =  0-*«/=l,M 


Fig.  1.  Magnetically  compensated  supersonic  beam. 


0146-9592/89/231324-03$2.00/0  ©  1989  Optical  Society  of  America 


December  1.  1989  /  Vol.  ll.  No.  ^3  /  OFI  iCS  LEITERS 


=  1  transition  at  690  nm  is  excited.  Without  magnetic 
compensation,  the  Doppler-broadened  linewidth  for 
this  configuration  is  found  to  be  ~480  MHz. 

To  reduce  the  Doppler  linewidth,  we  exploit  the 
low-temperature  narrow  velocity  distribution  of  the 
supersonic  beam. '  In  this  case,  the  Doppler  frequen¬ 
cy  shift  of  a  given  atom  is  approximately  wo^A,  where 
uo  is  the  average  stream  speed,  6  is  the  angle  between 
the  atom  direction  of  motion  and  the  perpendicular  to 
the  laser  beam,  and  A  is  the  optical  wavelength.  Since 
0  is  approximately  linear  in  distance  along  the  laser 
beam  axis,  the  Doppler  frequency  shift  is  also  approxi¬ 
mately  linear  with  position  along  the  laser  axis. 
Hence,  a  spatially  varying  2ieeman  shift  can  be  used  to 
tune  each  atom  into  resonance  with  the  local  Doppler- 
shifted  laser  frequency,  so  that  the  entire  interaction 
length  of  the  atomic  beam  can  be  made  to  resonate  at 
the  same  frequency.  The  field  produced  by  a  pair  of 
opposing  magnet  coils  (Fig.  1)  is  used  to  Zeeman  tune 
the  M  =  1  excited  state,  for  which  the  g  factor  is  3/2. 
The  aspect  ratio  for  these  coils  is  chosen  for  maximum 
field  gradient,  and  the  coils  are  water  cooled.  In  a 
typical  experiment  the  gradient  required  to  compen¬ 
sate  the  Doppler  broadening  is  only  approximately  50 
G/cm  for  the  Zeeman  tuning  rate  of  2  MHz/G. 

Figure  2  shows  the  absorption  enhancement  and 
line  narrowing  obtained  using  magnetic  compensation 
compared  with  that  obtained  without  compensation. 
In  this  case,  the  maximum  linear  absorption  is  limited 
to  50%  in  order  to  investigate  the  line  shape.  Absorp¬ 
tions  up  to  ~80%  were  obtained  in  this  system  at 
higher  oven  temperatures.  With  the  present  setup, 
nozzle  clogging  often  occurred  at  still  higher  oven  tem¬ 
peratures  when  Sr  was  used  as  the  active  species,  al¬ 
though  this  did  not  happen  for  Yb. 

In  order  to  understand  the  line  shape  for  the  com¬ 
pensated  beam,  the  path-averaged  linear  absorption 
coefficient  was  calculated  in  first-order  perturbation 
theory.  The  result  is  easily  understood  using  heuris¬ 
tic  arguments  and  is  given  by 


where  A  is  the  laser  detuning  and  the  angle-dependent 
linewidth  is  given  by 
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The  laser  beam  centered  at  a  distance  from  the 
source,  yj,  has  an  intensity  1/e  radius  R  and  a  wave 
vector  q.  The  supersonic  atomic  speed  distribution  is 
assumed  to  be  a  Gaussian  of  1/e  width  u  relative  to  the 
mean  speed  uo-  It  is  assumed  also  that  the  atomic 
density  does  not  vary  across  the  interaction  path 
length  for  the  small  angles  of  interest  here.  This  is 
consistent  with  the  flatness  of  the  experimental  un¬ 


compensated  absorption  line  shape  that  effectively 
measures  the  atomic  spatial  distribution.  The  Zee- 
man  tuning  rate  per  centimeter  is  d  =  {gtto/h){dBo/dz). 

Gaussian  velocity  and  starting  position  distribu¬ 
tions  are  assumed,  so  that  the  line  shape  in  Eq.  (2)  is  a 
sum  of  Gaussian  line  shapes  for  each  angle  6.  The 
spontaneous  linewidth  is  considered  negligible  com¬ 
pared  to  other  sources  of  line  broadening.  With  the 
above  definitions,  primary  contributions  to  the 
linewidth  arise  from  the  laser  beam  transit  rate  vq/R 
and  from  the  variation  of  the  Zeeman  shift  due  to  the 
distribution  of  the  starting  positions  of  atoms  across 
the  oven  nozzle,  which  is  taken  to  be  a  Gaussian  of  1/e 
width  ro.  For  a  point  source  this  contribution  vanish¬ 
es.  Angle-dependent  contributions  to  the  linewidth 
arise  from  the  Zeeman  frequency  chirp  0R0  as  atoms 
move  across  the  laser  beam  and  from  the  longitudinal 
residual  Doppler  width  qu0,  which  results  from  the 
projection  of  the  width  of  the  sp>eed  distribution  onto 
the  laser  beam  propagation  direction.  Note  that  for  a 
given  speed  the  Doppler  frequency  shifts  correspond¬ 
ing  to  different  angles  0  are  included  in  the  analysis 
through  the  uniform  angular  distribution.  These 
Do  'pler  shifts  are  compensated  by  the  magnetic  field 
gradient.  The  frequency-integrated  absorption  coef¬ 
ficient  is  just  the  usual  two-level  atom  result 
no(4T2k^2/A)  =  (W2)7,pontffoptno.  where  itopt  =  (S/2t)\‘^ 
and  no  is  the  ^Sr  atomic  density  in  the  interaction 
region.  For  an  interaction  path  length  I,  the  measured 
weak  field  transmitted  power  fraction  will  be 
exp(-a(A)l]. 

To  compare  the  calculated  and  experimentally  ob¬ 
served  absorption  line  shapes,  input  parameters  were 
calculated  for  the  measured  oven  temperature  of 
~1100*C  using  standard  one-dimensional  adiabatic 
expansion  results.'-^  The  total  oven  density  is  esti¬ 
mated  from  the  vapor  pressure  curve’'  to  be  n  =  70.2  X 
10*®/cc.  The  corresponding  stream  speed  is  uo  =  1.3  X 
10^  cm/sec.  Assuming  a  collision  cross  section  of  60 
and  a  collision  effectiveness  of  0.25,  one  obtains  a 
terminal  Mach  number  of  6.2  and  a  velocity  spread  of 


Pig.  2.  Uncompensated  (lower  trace)  and  compensated 
(upper  trace)  absorption  line  shapes  in  Sr.  The  maximum 
compensated  absorption  is  50%.  The  larger  (smaller)  peak 
is  the  ""Sr  (“Sr)  isotope. 
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Pig.  3.  Theoretical  compensated  and  uncompensated  ab¬ 
sorption  coefficients  versus  the  detunin"  in  units  of  ao  = 


u  =  2.3  X  10^  cm/sec.^  These  numbers  yield  a  residual 
Doppler  shift  of  qu  =  2.1  X  10^  rad/iec  per  radian  and 
an  uncompensated  Doppler  line  broadening  of  qvo  = 
1.2  X  10’°  rad/sec  per  radian.  For  yi  =  14  cm  the 
optimum  Zeeman  shift  per  unit  length  is  /3  =  qvo/yi  = 
0.84  X  10®  rad/sec  per  centimeter.  With  a  laser  beam 
of  0.5-mm  intensity  1/e  radius,  the  frequency  chirp  per 
radian  is  0R  =  0.42  X  10®  rad/sec,  which  is  negligible 
compared  to  the  residual  Doppler  shift  per  radian. 
Since  the  spontaneous  lifetime  is  long,  the  primary 
contribution  to  the  angle-independent  linewidth 
arises  from  the  oven  aperture  and  the  compensating 
field  gradient.  For  a  nozzle  radius  of  ro  =  75  urn,  this 
linewidth  is  ffro  =  0.56  X  10^  rad/sec.  With  the  com¬ 
pensating  field  off,  the  corresponding  linewidth  is 
transit  time  limited  with  Vo/R  =  0.26  X  10^  rad/sec  for 
a  laser  beam  intensity  1/e  radius  of  0.5  mm.  Using  these 
parameters,  one  obtains  the  results  shown  in  Fig.  3. 

The  maximum  calculated  absorption  at  line  center 
can  be  used  to  estimate  the  atomic  density,  which  is  0.8 
X  10~*  (ir/2)7,p„nt<ropt»io-  For  the  measured  50%  ab¬ 
sorption,  one  can  estimate  the  atomic  density  using  al 
=  0.69, 1  =  4.4  cm,  <r„pi  =  2.3  X  10"®  cm®  at  690  nm,  and 
7»pont  =  4.46  X  10''  sec"*.  This  yields  no  =  1.2  X  10**/ 
cc.  For  comparison,  the  total  Sr  density  calculated 
from  the  oven  density  at  a  distance  L  =  14  cm  from  the 
nozzle  (assuming  that  the  atoms  travel  in  straight  lines 
after  the  adiabatic  expansion  ends)  is  0.157(2ro/L)®n  = 
1.27  X  10“.®  Using  the  ®®Sr  abundance  of  83%,  we 
have  no  =  1.05  X  10“. 

The  calculated  ratio  of  peak  compensated  to  un¬ 
compensated  absorption  is  29.0:1,  which  compares 
well  with  the  measured  value  of  29:4:1.  Note  that  in 
the  Doppler-broadened  limit,  collimation  of  an  un¬ 
compensated  beam  will  reduce  the  linewidth  without 
changing  the  uncompensated  absorption.  Hence,  the 
improvement  over  a  collimated  atomic  beam  with  the 
same  source  consumption  and  linewidth  is  roughly  of  a 
factor  of  29.  Higher  absorption  coefficients  and  nar¬ 
rower  linewidths  can  be  obtained  at  a  higher  source 
density.  However,  collisions  outside  the  nozzle  ulti¬ 


mately  increase  the  effective  source  size  and  thereby 
increase  the  linewidth. 

While  the  calculated  compensation  ratio  and  densi¬ 
ty  agree  well  .«ith  the  experimental  values,  the  line 
shape  obtained  using  the  parameters  based  on  the 
oven  temperature  and  oven  diameter  is  qualitatively 
different  from  that  measured.  Better  agreement  is 
obtained  by  enlarging  the  effective  nozzle  diameter  by 
a  factor  of  6  to  take  into  account  collisions  in  the 
expansion  region,  thus  increasing  the  angle-indepen¬ 
dent  linewidth.  The  residual  Doppler  width  is  then 
reduced  by  assuming  a  higher  Mach  number  of  12.4  (a 
larger  collision  cross  section) .  This  decreases  the  ped¬ 
estal  on  the  calculated  line  shape,  and  the  HWHM  can 
be  made  to  agree  with  the  experimentally  observed 
value  of  approximately  18  MHz.  Since  parameters 
are  chosen  to  obtain  agreement  with  the  compensated 
and  uncompensated  maximum  absorption  at  line  cen¬ 
ter,  the  calculated  atomic  density  obtained  in  the  in¬ 
teraction  region  remains  as  stated  above. 

The  compensated  beam  system  demonr+rated  here 
provides  a  reasonably  general  technique  for  obtaining 
high  atomic  densities  and  long  interaction  path 
lengths  for  nonlinear  optics  and  spectroscopic  experi¬ 
ments  while  a  narrow  linewidth  is  maintained.  It  can 
be  applied  whenever  tbe  Doppler  shift  and  the  posi¬ 
tion  are  correlated.  In  addition  to  supersonic  beams, 
any  type  of  accelerated  beam,  for  which  the  fractional 
velocity  spread  is  small,  can  be  compensated  to  reduce 
the  spectral  linewidth.  We  have  applied  the  compen¬ 
sated  beam  technique  to  the  *So  -*•  transition  of 
*®^Yb,  where  the  sample  was  easily  made  optically 
thick  although  the  spontaneous  lifetime  is  875  nsec. 
While  this  technique  is  effective  for  only  one  M  level 
and  propagation  direction,  it  may  have  important  ap¬ 
plications  to  atom-cavity  experiments  requiring  a  high 
atomic  density  of  two-level  atoms  and  to  nonlinear- 
optics  experiments  in  general.  In  addition,  for  weak- 
line  spectroscopy,  this  method  may  provide  a  means  of 
obtaining  substantial  signal  enhancement. 

This  research  is  supported  by  the  Rome  Air  Devel¬ 
opment  Center  under  contract  Fl9628-88-K-()012. 
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We  report  the  observation  of  continuous  spatial  photon  echoes  in  a  magnetically  compensated  supersonic  beam  of 
atomic  samarium.  Atoms  cross  two  spatially  separated  cw  laser  pump  Helds  to  produce  an  intense  cw  echo 
downstream.  In  contrast  to  previous  studies  of  continuous  coherent  radiation  in  separated  fields,  the  signals  are 
readily  observable  with  the  unaided  eye.  The  echo  signal,  which  is  generated  in  a  magnetic-Held  gradient,  differs 
from  usual  rephasing  phenomena  owing  to  a  quadratic  time  dependence  for  the  optical  phase  in  the  frame  of  the 
moving  atoms. 


Since  the  first  observation  of  the  photon  echo,*  optical 
rephasing  techniques  have  found  a  variety  of  applica¬ 
tions  in  studies  of  relaxation  phenomena  in  vapors  and 
condensed-matter  systems.  Recently,  rephasing 
techniques  have  been  applied  to  optical  processing'^ 
and  information  storage.  *  '’  Rephasing  techniques 
also  have  found  important  application  in  developing 
optical  Ramsey  resonance  techniques  in  atomic 
beams.® 

In  this  Letter  we  report  the  observation  of  a  new 
type  of  continuous  spatial  photon  echo  in  a  magneti¬ 
cally  compensated  supersonic  beam^  of  samarium  at¬ 
oms.  By  using  a  magnetic-field  gradient  to  compen¬ 
sate  the  Doppler  shifts  of  the  diverging  supersonic 
beam,  a  high  product  of  density  and  path  length  (nl)  is 
obtained.  In  this  case,  substantial  optical  coherence 
is  spatially  transferred  in  the  experiments,  leading  to 
large  coherent  signals  that  are  radiated  and  directly 
detected  downstream  from  two  pump  regions  that  the 
atoms  traverse.  Continuous  coherent  radiation  by  at¬ 
oms  that  have  traversed  spatially  separated  pump 
fields  leads  to  signals  proportional  to  the  square  of  the 
number  of  atoms  and  so  is  closely  related  to  superra¬ 
diance  in  extended  media.®  Coherent  radiation  in 
separated  fields  has  been  studied  previously  in  vapors® 
without  the  use  of  a  magnetically  compensated  beam. 
However,  heterodyne  detection  was  needed  to  mea¬ 
sure  the  10" ’  ’’-W  signals  that  were  obtained. 

In  our  experiments  (Fig.  1)  a  supersonic  atomic 
beam  traverses  two  continuous  pump  regions.  These 
regions  play  the  roles  of  the  usual  )r/2  and  ir  pulses  in 
temporal  rephasing  phenomena  as  seen  in  the  atom 
frame.  Hence,  inhomogeneous  dephasing  of  the  mac¬ 
roscopic  polarization  during  the  traversal  time  be¬ 
tween  the  two  pump  beams  is  reversed  after  the  sec¬ 
ond  pump  field,  leading  to  a  large  macroscopic  polar¬ 
ization  downstream  from  the  second  pump  region. 
This  idea  has  been  exploited  previously  in  three-beam 
optical  Ramsey  resonance  experiments  to  eliminate 
optical  dephasing.®  Unlike  that  in  previous  experi¬ 
ments,  which  measure  a  small  absorption  or  fluores¬ 
cence,  our  system  employs  a  magnetically  compensat¬ 
ed  supersonic  beam  to  achieve  a  long  interaction  path 
length  in  the  medium  with  a  narrow  Doppler  band¬ 


width.  This  permits  a  large  optical  polarization  to  be 
obtained  so  that  the  echo  Held  generated  by  the  atomic 
medium  0.5  mm  downstream  from  the  pump  region  is 
visible  with  the  unaided  eye,  in  contrast  to  previous 
experiments  on  continuous  radiation  in  spatially  sepa¬ 
rated  fields.® 

The  system  employs  a  supersonic  samarium  beam 
that  emerges  from  a  IbO-^m  pinhole  drilled  in  the  side 
of  an  electrically  heated  tantrum  tube.  This  beam  is 
collimated  to  a  1-cm  aperture,  where  it  crosses  two  cw 
linearly  polarized  pump  fields  of  powers  of  a  few  hun¬ 
dred  microwatts.  Using  a  SO-tim  waist  and  a  2-mm 
height  for  the  pump  beams,  we  achieve  a  transit  time 
bandwidth  of  5  MHz.  The  ^Fo  -*  ’D  i  transition  at  599 
nm  is  excited  for  which  the  radiative  lifetime  is  1.7 
lisec.  Without  magnetic  compensation,  the  Doppler 
width  is  ~1()0  MHz,  and  few  atoms  are  excited  by  the 
pump  fields.  To  reduce  the  Doppler  width,  the  Dopp¬ 
ler  shifts  are  compensated  by  means  of  two  opposing 
magnet  coils,  which  provide  a  Zeeman  shift  that  varies 
linearly  with  position  along  the  laser  beam  axis.^  The 


Fig.  1 .  Magnetically  compensated  beam  showing  the  pump 
beams  and  the  echo. 
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Zeeman  shift  gradient  is  chosen  so  that  the  'Di,  Af  =  1 
state  (g  factor  2.83)  is  resonant  with  the  Doppler- 
shifted  laser  frequency  at  each  position  along  the  laser 
beam  axis.  This  method  is  closely  related  to  that 
employed  previously  for  laser  cooling,  in  which  a  spa¬ 
tially  varying  Zeeman  shift  compensates  for  the  Dopp¬ 
ler  detuning  that  arises  as  the  atoms  are  slowed. For 
optimum  magnetic-field  gradient,  the  Doppler  width 
is  reduced  to  5  MHz,  so  that  the  transit  time  band¬ 
width  is  adequate  to  excite  most  atoms  and  a  large 
echo  signal  is  obtained.  By  contrast,  without  magnet¬ 
ic  compensation,  echo  signals  are  reduced  by  more 
than  two  orders  of  magnitude.  The  signals  are  ob¬ 
served  through  a  crossed  linear  polarizer  so  that  the 
circularly  polarized  echo  and  the  circularly  polarized 
free-induction  decay  signals  induced  by  the  pump 
fields  are  of  similar  intensity,  while  the  pump  fields 
are  blocked.  These  signals  are  imaged  onto  a  Reticon 
diode  array  (spacing  25  /im).  The  laser  is  servo  stabi¬ 
lized  to  the  atomic  resonance  frequency  to  eliminate 
long-term  drift.  A  typical  echo  signal  is  shown  in  Fig. 
2.  Note  that  the  first  rf2  pulse  is  more  effective  than 
the  second  ir  pulse  in  producing  free-induction  decay 
signals,  so  that  the  lower-intensity  pulse  yields  the 
larger  transmitted  signal. 

The  presence  of  a  magnetic- field  gradient  in  addi¬ 
tion  to  residual  Doppler  broadening  leads  to  a  qua¬ 
dratic  time  dependence  for  the  phase  of  the  macro¬ 
scopic  polarization  in  the  atom  frame.  This  alters  the 
variation  of  the  echo  intensity  with  position  along  the 
atomic  propagation  axis  from  that  normally  expected 
for  a  two-pulse  echo,  as  shown  in  Fig.  3.  For  narrow 
pump  beams,  the  basic  features  of  echo  formation  in  a 
compensated  beam  can  be  derived  by  using  heuristic 
arguments  as  follows.  We  take  the  atomic  beam  axis 
in  the  y  direction,  with  the  center  of  the  first  pump 
region  at  y  =  yi  and  that  of  the  second  pump  region  at 
y  =  y2-  For  an  atom  of  velocity  v,  optical  polarization 
is  created  in  the  first  interaction  region  at  time  t\  =  yj 
Uy  assuming  that  the  atom  is  at  the  oven  nozzle  at  t  =  0. 
The  polarization  acquires  a  phase  that  is  just  the  time 
integral  of  the  instantaneous  optical  frequency.  The 
instantaneous  optical  frequency  differs  from  the  natu- 
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Fig.  2.  Continuous  spatial  photon  echo.  The  pump  beams 
viewed  through  a  crossed  linear  polarizer  appear  on  the  left. 
The  echo  signal  is  the  rightmost  peak. 
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Fig.  3.  Decay  of  the  spatial  echo  intensity  with  position. 
The  data  points  are  obtained  by  translating  the  position  of 
the  IT  beam  relative  to  the  Tt/2  beam. 


ral  frequency  owing  to  the  Doppler  shift  kv^  and  the 
Zeeman  shift  —az,  where  the  z  axis  is  taken  as  the  laser 
beam  axis  and  k  =  2jr/X  is  the  optical  wave  vector. 
The  Zeeman  shift  —az  =  —(giJ.o/h){dB/dz)z  is  assumed 
to  vary  linearly  with  position  along  the  z  axis.  The 
atomic  position  at  time  f  is  y  =  Pyt  and  z  =  v^t.  Hence, 
the  frequency  of  the  polarization  for  times  f  >  1 1  is  A(t) 
=  kVi  —  aVit.  The  net  phase  accumulated  for  an  atom 
traveling  between  the  two  regions  is  then  vr2i  = 
S',\  dtA(t)  =  kVt(t2  —  ti)  -  If  we 

assume  that  in  traversing  the  second  region  the  atomic 
Bloch  vector  is  rotated  through  a  ir  angle,  the  accumu¬ 
lated  phase  is  reversed.  Following  this,  the  atom 
propagates  downstream  to  a  position  y,  where  its  accu¬ 
mulated  phase  is  v’Cy)  =  /[,  dtA(()  —  <P2\,  with  t  =  y/Vy. 
Hence,  the  phase  of  the  optical  polarization  at  position 
y  is 


<p(y)  =  k^(y  -2y2  +  yi)  -  ~^{y^  -  2y^  -I-  y,^). 

(1) 

For  the  supersonic  beam,  the  velocity  component  c,  is 
distributed  as  a  Gaussian  with  a  narrow  spread  u 
about  the  mean  speed  on-  The  net  optical  dipole  mo¬ 
ment  is  obtained  by  integrating  the  atomic  dipole  mo¬ 
ment  over  the  speed  distribution  and  over  the  interac¬ 
tion  volume.  The  path  length  along  the  laser  propa¬ 
gation  drection  is  proportional  to  the  angle  6  =  vjvy. 
The  dipole  moment  then  is  proportional  to 

/  f  /  MV  ■  .L  ^  f""’"  do 

(exp[np(y)J>  =  sinfli, - - -  — — 

where  the  pulse  areas  in  the  atom  frame  are  (S>\  and  <t>2 
for  pump  regions  1  and  2,  respectively.  In  writing  Eq. 
(2),  we  assume  that  the  residual  Doppler  detuning  is 
small  compared  with  the  transit  time  bandwidth  and 
the  Rabi  frequencies.  The  integral  over  0  is  readily 
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carried  out.  For  simplicity,  we  assume  that  the  mag¬ 
netic-field  gradient  is  chosen  for  ideal  frequency  com¬ 
pensation  in  the  first  pump  region  aty  =  yi,  so  thatd<p/ 
dy  =  0aty  =  yi.  Then  «  = /fOii/yi.  Since  the  pump 
beams  and  echo  occur  within  approximately  a  milli¬ 
meter  of  yi,  withy]  =  14  cm,  wedefiney  =yi  -l-fandy2 
=  yi  +  (o.  In  this  case,  for  7r/2  and  ir  pulses  in  pump 
regions  1  and  2,  the  dipole  moment  is  then  proportion¬ 
al  to 

,  ,  f*  A/nd.t  „ 

(exp[iV(y)])  = - exp(-Af„'x-) 

}-<x.  ^ 


U  -  2<  ,) 


('  -  2(o'  1 


X  +  1 


2y,  X  -I-  1 


<7^n 


U  -  2f,) 


r  ■  -  2<2‘  1 


“  X  -h  1  2y|  X  +  1 


(3) 


where  x  =  {l\  —  upl/on  and  A/p  =  vn/u  is  the  Mach 
number.  From  Eq.  (3),  it  is  clear  that  the  spread  in 
the  speed  distribution  leads  to  dephasing,  and  rephas¬ 
ing  tends  to  occur  at  position  <  =  2(2  (i.e.,  at  twice  the 
separation  between  the  two  pump  regions).  However, 
the  quadratic  dependence  of  the  phase  on  distance 
that  arises  from  the  magnetic-field  gradient  favors 
rephasing  at  a  point  where  =  2<2^.  Indeed,  as  u  -►  0 
at  high  Mach  numbers,  x  -*  0,  and  the  rephasing  point 
is  (  =  V  2<2-  For  pump  region  separations  (2  of  ~1  mm, 
the  quadratic  position  dependence  of  phase  of  the  sine 
function  dramatically  alters  the  dependence  of  the 
echo  intensity  on  position  from  the  first  pump  region 
from  that  expected  for  usual  two-pulse  echoes.  Our 
initial  measurements  of  echo  intensity  versus  pump- 
beam  separation  are  in  qualitative  agreement  with  Eq. 
(3)  assuming  a  Mach  number  Mo  =  5,  after  averaging 
over  a  spread  of  <2  in  order  to  account  for  the  finite 
pump-beam  diameter  as  shown  in  Fig.  3.  However, 
complete  quantitative  agreement  has  not  yet  been 
obtained. 

In  conclusion,  we  have  demonstrated  a  method  for 
preparing  and  transporting  a  large  macroscopic  opti¬ 
cal  coherence  over  millimeter  distances,  so  that  optical 
coherence  can  be  observed  in  a  region  free  of  all  pump 
fields.  While  the  local  frequency  compensation  tech¬ 
nique  demonstrated  in  Ref.  7  can  be  used  to  obtain 
large  signals,  the  observed  decay  of  the  echo  intensity 
with  increasing  distance  shows  that  the  quadratic  time 
dependence  of  the  optical  phase  (for  the  present  mag¬ 
netic-field  gradient  configuration)  limits  the  maxi¬ 
mum  distances  over  which  coherence  can  be  trans¬ 
ported  without  dephasing. 


The  method  demonstrated  in  this  Letter  is  ideally 
suited  to  developing  spatial  field  convolution  and  cor¬ 
relation  techniques.  Such  techniques  have  been  ex¬ 
plored  previously  using  four-wave  mixing,'*  except 
that  the  individual  interactions  with  the  input  fields 
were  not  spatially  separated.  As  applied  to  optical 
processing,  our  techniques  are  the  spatial  analogs  of 
the  methods  used  in  the  time  domain  based  on  stimu¬ 
lated  echoes.^'*^  The  present  experiments  achieve  a 
spatial  resolution  of  25  >im  limited  by  the  Reticon 
diode  array.  Much  higher  spatial  resolution  (up  to 
the  diffraction  limit)  is  attainable  by  magnification  of 
the  echo  image  and  increasing  the  bandwidth  of  the 
medium,  which  is  fully  adjustable.  In  addition,  opti¬ 
cal  Bloch  vector  noise  can  be  studied  as  a  function  of 
the  Bloch  angle'^  by  using  a  cw  pump  fields  and  mea¬ 
suring  the  phase-dependent  noise  in  the  cw  echo  sig¬ 
nals.  Finally,  the  method  is  well  suited  to  study  the 
evolution  of  optical  coherence  in  cavities''*  in  which 
the  emitted  field  can  be  measured  directly  and  contin¬ 
uously  by  optical  methods. 

This  research  is  supported  by  the  Rome  Air  Devel¬ 
opment  Center  under  contract  F19628-88-K-0012. 
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A  method  for  measuring  atomic  spatial  distributions  is  suggested,  based  on  Raman  transitions  in  an  applied  force 
Held.  This  technique  eliminates  the  Doppler  effect  and  provides  high-frequency  resolution  for  long-lived  states, 
leading  to  position  resolution  that  is  fundamentally  limited  by  the  uncertainty  principle.  The  effects  of  atomic 
motion  are  discussed,  and  a  simplified  analysis  of  the  method  is  presented. 


Recently  there  has  been  substantial  progress  and  in¬ 
terest  in  the  creation  of  atomic  distributions  that  vary 
over  small  length  scales  by  diffraction,*  periodic  spa¬ 
tial  modulation, 2  channeling,’’  focusing, and  cooling.® 
In  general,  such  distributions  exhibit  momentum- 
space  coherence  and  are  of  both  practical  and  funda¬ 
mental  interest.  These  developments  have  been 

urred  on,  in  part,  by  the  advent  of  practical  laser 
^ooling  techniques.®  Important  applications  include 
atomic  interferometry  and  the  creation  of  submicro¬ 
meter  structures  by  atomic  deposition.  While  tech¬ 
niques  exist  to  create  suboptical-wavelength  atomic 
spatial  modulation  and  interference,  with  the  excep¬ 
tion  of  the  technique  described  in  Ref.  3,  which  is 
described  below,  methods  for  detection  have  been  lim¬ 
ited  principally  to  hot  wires,  so  far  as  the  author  is 
aware. 

In  this  Letter  a  class  of  techniques  is  suggested  that 
may  prove  useful  in  achieving  much  higher  spatial 
resolution  of  moving  atoms  by  using  optical  fields  and 
that  examine  interesting  features  of  quantum-me¬ 
chanical  position  measurement. 

The  method  for  making  the  measurements  is  in 
principle  quite  simple.  Suppose  that  an  atom  in  a 
state  (0;  Fig.  1)  enters  a  region  of  potential  that,  for 
simplicity,  varies  linearly  in  one  spatial  direction  x  and 
affects  only  level  2.  Then  the  potential  can  be  written 
as 

V(X)  =  -(X  -  XgfF,  (1) 

where  f  is  a  constant  force.  A  number  of  techniques 
can  be  used  to  generate  such  a  potential,  some  of  which 
are  described  below.  The  basic  idea  is  to  use  the  fact 
that  for  a  transition  from  level  0  to  2,  the  spatially 
varying  potential  correlates  the  resonance  frequency 
of  the  atom  with  its  position,  so  that  optical  fields  of 
sufHcient  frequency  resolution  can  select  out  atoms  at 
one  particular  position,  analogous  to  a  nuclear  mag¬ 
netic  resonance  imaging  technique.  Doppler  frequen¬ 
cy  shifts,  which  broaden  the  line  for  one-photon  tran¬ 
sitions,  are  readily  eliminated  by  using  a  copropagat¬ 
ing  wave  Raman-induced  optical  resonance^  (Fig.  1)  or 
any  other  Doppler-free  coherent  process  connecting 
long-lived  states.  In  the  Raman  transition,  it  is  the 


difference  between  the  frequencies  of  two  optical 
flelds  that  must  resonate  with  the  atomic  system,  and 
the  Doppler  shift  approximately  cancels,  so  the  posi¬ 
tion  measurement  can  be  made  independent  of  veloci¬ 
ty.  The  above  potential  shifts  the  final  atomic  state  2 
and  tunes  the  Raman  resonance  frequency  to  make  it 
spatially  dependent  with  a  frequency  gradient  F/h  = 
dco/dx.  In  practice,  the  optical  field  difference  fre¬ 
quency  can  be  taken  as  the  unperturbed  initial-to- 
final  spacing,  so  that  the  above  potential  selects  out 
atoms  at  position  x  =  xq,  the  position  at  which  the  net 
frequency  shift  of  the  spatially  varying  potential  is 
zero.  This  position  is  readily  moved  by  varying,  for 
example,  a  small  uniform  magnetic  field,  which  is  su¬ 
perimposed  upon  the  spatially  varying  potential. 

An  importwt  feature  of  the  Raman  technique  is 
that  it  corresponds  to  the  absorption  of  a  photon  from 
one  optical  field  and  the  emission  of  a  photon  into 
another,  copropagating,  field.  Hence  the  net  momen¬ 
tum  imparted  by  the  Raman  optical  fields  to  the  atom 
is  approximately  zero,  since  the  difference  frequencies 
of  interest  here  are  small.  This  has  the  interesting 
consequence  that  the  accuracy  of  the  measurement  is 
fundamentally  limited  by  the  Heisenberg  uncertainty 
principle,  even  when  the  Raman  fields  propagate 
along  the  measurement  axis.  Before  we  discuss  other 
limiting  factors  that  are  due  to  the  atomic  motion,  it  is 
useful  to  consider  this  limit.  First,  consider  the  fre¬ 
quency  resolution  required  to  obtain  a  spatial  resolu- 


Fig.  1.  Three-level  system  and  Raman  fields.  Atoms  ini¬ 
tially  are  in  state  0. 
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tion  Ax.  The  resonance  frequency  shift  across  the 
spatial  resolution  length  scale  must  equal  the  frequen¬ 
cy  resolution: 

FLx/h  =  1/T.  (2) 

Since  the  initial  and  final  atomic  states  are  taken  as 
long-lived  excited  states  or  ground  states,  the  frequen¬ 
cy  resolution  will  be  limited  by  the  interaction  time  T 
of  the  atom  with  the  optical  field.  Hence  Ax  =  h/ 
(FT).  Raman  transitions  can  occur  at  any  time  in  the 
Raman  region,  so  the  maximum  momentum  imparted 
by  the  potential  is  just  Pni«i  =  FT.  Estimating  Ap  ~ 
Pmax  yields  AxAp  =  h,  which  determines  the  funda¬ 
mental  limit.  Thus  the  resolution  accuracy  is  deter¬ 
mined  by  the  spatially  varying  potential  according  to 
the  uncertainty  principle,  and  no  excess  momentum  is 
imparted  by  the  measurement  light  field. 

As  described  below,  a  number  of  techniques  are 
available  for  producing  spatially  varying  level  shifts  of 
10*°  Hz/cm  over  millimeter  distances,  with  much  larg¬ 
er  shift  gradients  obtainable  for  smaller  maximum 
distances.  Assuming  a  transit-time  limited  resolution 
of  100  kHz,  the  corresponding  spatial  resolution  is  0.1 
pm.  By  using  a  Raman-induced  optical  resonance 
technique,  it  is  easy  to  obtain  much  better  frequency 
resolution  than  100  kHz  by  increasing  T.  However,  as 
T  increases,  the  atomic  motion  limits  the  spatial  reso¬ 
lution.  In  the  simplest  case,  the  spatial  resolution  will 
be  limited  by  the  velocity  of  the  atoms  along  the  mea¬ 
surement  axis.  To  determine  the  natural  scale  of  res¬ 
olution,  we  equate  the  interaction  time  T  with  the 
transit  time  of  the  atoms  across  the  resolution  length. 
If  we  take  the  atomic  beam  to  be  principally  propagat¬ 
ing  along  the  y  axis  with  velocity  Uy  and  with  a  compo¬ 
nent  of  velocity  along  the  measurement  x  axis,  u,,  then 
the  time  to  cross  the  resolution  length  is  T  -  Ax/vx. 
The  corresponding  frequency  resolution  is  1/T  =  o*/ 
Ax.  Using  this  in  Eq.  (2)  and  solving  for  Ax  yields 

Axvel  =  (huJF)'^''^.  (3) 

A  frequency  gradient  of  10'°  Hz/cm  gives  Axvf.l  =  0.04 
tira[vx  (cm/sec)]  For  a  supersonic  beam  of  velocity 
Vy  =  10®  cm/sec,  modest  collimation  of  1  mrad  would 
yield  Cj  =  100  cm/sec  and  a  spatial  resolution  of  M).4 
pm.  The  optical  laser  beam  width  d  must  then  be 
chosen  for  consistency  such  that  the  laser  beam  transit 
time  equals  the  transit  time  across  the  spatial  resolu¬ 
tion  length:  d  =  Ax(vy/ux)  =  0.4  pm/(l  mrad)  =  0.4 
mm.  In  this  way,  the  longest  transit  time  possible  is 
used,  for  which  the  atom  does  not  travel  any  farther 
across  the  frequency  gradient  than  the  resolution 
length.  For  this  reason,  optical  lasers  are  particularly 
well  suited  to  limiting  the  interaction  time  by  focusing. 

If  the  velocity  along  the  measurement  axis  is  suffi¬ 
ciently  small,  the  acceleration  due  to  the  spatially 
varying  potential  is  a  limiting  factor.  In  this  case, 
equating  the  resolution  length  with  the  position 
change  due  to  the  acceleration  yields  the  maximum 
interaction  time  T.  Inserting  this  into  Eq.  (2)  yields 
the  optimum  resolution  as 

A'tACCF.i.  =  (hV2MF)^'\  (4) 

For  consistency,  the  laser  beam  width  is  chosen  as  d  = 


VyT  to  limit  the  interaction  time.  In  this  case,  the 
momentum  imparted  by  the  force  limits  the  measure¬ 
ment  accuracy,  so  this  limit  is  in  some  sense  the  uncer¬ 
tainty  limit.  For  a  sodium  atom,  a  10*°-Hz/cm  fre¬ 
quency  gradient  yields  a  resolution  of  0.06  pm.  For  a 
velocity  perpendicular  to  the  laser  of  10®  cm/sec,  the 
corresponding  laser  beam  width  is  2.5  mm.  The  con¬ 
dition  for  the  velocity  along  the  measurement  axis  to 
dominate  the  resolution  is  for  Axvel  —  Axaccel-  This 
requires  that  the  velocity  o,  be  such  that  Uj  i 
Ax^ccel^/*’  which  is  equivalent  to  Vx  ^ 
/i/(2MAxaccel)-  For  the  above  conditions  this  yields 
Ox  >  2  cm/sec. 

Omitted  in  the  above  discussion  are  the  limits  im¬ 
posed  by  the  atomic  distribution  itself.  In  general, 
any  atomic  beam  will  have  some  finite  transverse  ve¬ 
locity  spread.  In  a  state  possessing  momentum-space 
coherence,  the  different  kinetic  energies  correspond¬ 
ing  to  different  momentum  components  will  cause  pe¬ 
riodic  spatial  modulation  to  wash  out  when  the  (inco¬ 
herent)  transverse  velocity  averaging  is  performed. 
This  problem  has  been  solved  by  having  the  atoms 
traverse  multiple  separated  laser  fields  to  initiate  re- 
phasing.^  If  the  atomic  pattern  varies  more  slowly 
than  Axvel.  the  rephasing  length  will  be  longer  than 
the  laser  beam  width.  A  more  fundamental  limit  is 
due  to  wave-pocket  spreading.  In  this  case,  an  atomic 
distribution  that  varies  over  a  length  a  will  spread  a 
distance  of  order  a  over  a  time  Tmax  such  that  (h/ 
2Ma^)  Tmax  ~  1.  The  transit  time  across  the  resolu¬ 
tion  length  must  be  limited  to  T max-  Using  this  in  Eq. 
(2)  leads  to  a  resolution  limit  of  Axspread  = 
(Axaccel)^/o^-  Hence,  as  long  as  the  scale  of  spatial 
variation  is  larger  than  Axaccel.  wave-packet  spread¬ 
ing  during  the  measurement  poses  no  problem. 

Spatially  varying  frequency  shifts  can  be  produced 
by  a  variety  of  methods,  the  most  obvious  of  which  is 
the  use  of  Stark  or  Zeeman  shifts.  The  latter  have 
been  explored  extensively  for  atomic  beam  deflection 
applications.®  Field  gradients  larger  than  1(1*  G/cm 
are  obtained  by  using  hollow  water-cooled  two-wire 
configurations  and  millimeter  wire  spacings,  with  at¬ 
oms  propagating  parallel  to  the  wires.  In  order  to 
obtain  a  level  shift  that  varies  along  one  direction 
without  violating  the  divergence  or  curl  conditions,  a 
uniform  bias  field  is  applied,  either  as  part  of  the  two- 
wire  field  or  as  an  additional  field.  The  effect  of  this 
uniform  field  is  to  suppress  the  component  of  the 
spatially  varying  field  perpendicular  to  the  bias  field 
so  that  its  effect  on  the  level  shift  is  of  second  order  in 
distance.  For  sufficiently  large  bias  fields,  the  level 
shift  due  to  the  unwanted  field  variation  over  the  larg¬ 
est  distances  measured  can  be  made  smaller  than  the 
linear  variation  in  the  level  shift  over  the  resolution 
length. 

An  interesting  alternative  for  obtaining  spatially 
varying  level  shfits  is  to  exploit  the  light  shift  as  shown 
in  Fig.  2.  This  has  been  demonstrated  for  the  case  of 
atoms  channeled  in  a  nonresonant  standing  wave. 
Spatial  variation  of  the  light  shift  due  to  the  trapping 
light  field  was  used  to  deduce  the  periodic  variation  of 
atomic  density  from  the  one-photon  absorption  spec¬ 
tra.®  For  position  measurement  (Fig.  2),  two  nonreso¬ 
nant  linearly  polarized  laser  beams  of  wave  vector  q. 
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Fig.  2.  Spatial-frequency  gradient  based  on  the  light  shift 
of  two  overlapping  laser  beams.  Atoms  propagate  into  the 
page  along  the  y  axis.  Hi  and  02  are  the  Raman  laser  fields. 


intersect  at  an  angle  2d  and  are  used  to  shift  the  final 
state  2.  The  light  shift,  which  is  proportional  to  the 
intensity  in  lowest  order,  has  a  spatial  variation  <x  1  -h 
cos(29,x  sin  d).  Physically,  a  force  results  from  inter¬ 
actions  in  which  an  off-resonant  state  (state  2  in  this 
case)  absorbs  momentum  from  one  beam  and  emits 
into  the  other,  so  that  the  net  momentum  transferred 
to  the  atom  during  the  interaction  time  T  is  just  FT,  as 
in  the  above  discussion.  By  adjusting  the  intersection 
angle,  the  length  scale  over  which  the  light  shift  varies 
can  be  conveniently  controlled.  For  reasonable  de¬ 
tunings  and  laser  powers,  very  large  frequency  gradi¬ 
ents  are  obtainable  over  small  maximum  spatial  di¬ 
mensions.  For  example,  a  maximum  light  shift  of 
only  10  MHz  over  lO-^m  distances  would  3deld  fre¬ 
quency  gradients  of  10*®  Hz/cm.  For  measurement  of 
periodic  spatial  patterns  that  vary  at  multiples  of  the 
optical  wave  vector,  a  standing-wave  light  shift  can  be 
used  to  obtain  gradients  of  10*^  Hz/cm  at  moderate 
laser  intensity.  The  resonance  position  a;o  can  be  var¬ 
ied  near  the  maximum  derivative  point  either  by 
means  of  a  superposed  uniform  Zeeman  field  or  by 
tuning  the  Raman  laser  field  difference  frequency. 
Finally,  the  system  can  be  calibrated  by  using  atomic 
distributions  of  known  periodic  spatial  variation  pro¬ 
duced  by  any  convenient  method,*'®  and  the  shifting 
fields  can  be  stabilized  by  using  an  atomic  resonance. 
It  is  expected  that,  for  long-lived  atomic  states,  high 
spatial  resolution  can  be  obtained  in  this  way,  provid¬ 
ed  that  the  interaction  time  is  limited  to  the  optimum 
value,  as  discussed  above. 

Detailed  calculations  that  support  the  simple  argu¬ 
ments  discussed  above  have  been  carried  out  for  the 
level  scheme  of  Fig.  1,  treating  the  center-of-mass  mo¬ 


tion  for  the  atom  quantum  mechanically  and  including 
a  constant  force  that  acts  on  level  2  only.®  An  impor¬ 
tant  feature  is  that  the  detection  system  for  the  final 
state  need  not  have  any  spatial  resolution.  In  general, 
any  convenient  means  for  final-state  detection,  in¬ 
cluding  fluorescence,  photon  burst,  and  photoioniza¬ 
tion,  is  applicable.  This  permits  single-atom  detec¬ 
tion  schemes  and  either  ground  states  or  reasonably 
long-lived  excited  states  to  be  employed  for  the  mea¬ 
surements. 

The  maximum  state  2  count  rate  is  given  by  the 
product  of  the  Raman  transition  probability,  ijr,  the 
detection  efficiency,  i)d,  the  flux  in  state  0,  nvy,  and  the 
effective  interaction  area  that  is  the  product  of  the 
atomic  beam  width  perpendicular  to  the  measurement 
axis.  It,  and  the  spatial  resolution,  Axopt-  For  n  = 
10®/cm®  and  Vy  =  10®,  1*  =  0.1  mm,  Axopt  =  0.1  nm,  the 
count  rate  is  Tj^nj^/lOVsec.  Taking  ns  =  0.1  and  no  = 
0.1  yields  a  count  rate  of  10®/sec. 

Peuticularly  interesting  atomic  systems  for  mea¬ 
surements  are  the  ^F)  ’’Fo  transition  in  *®^Sm  and 
the  *So  transitions  of  *’^ Yb  and  ®®Sr.  The  latter 
two  have  a  long  radiative  lifetime  and  a  J  =  0  ground 
state  with  no  nuclear  spin.  They  can  be  used  by  ex¬ 
ploiting  the  excited-state  splittings  and  should  be  im¬ 
portant  for  atomic  interference  experiments  for  which 
unwanted  level  shifts  before  the  detection  region 
would  pose  a  problem. 
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A  simple  quantum  theory  of  recently  suggested  optical  techniques  for  ultrahigh-resolution  posi¬ 
tion  measurement  and  localization  of  moving  atoms  in  beams  is  presented.  Both  the  internal  and 
center-of-mass  motion  are  treated  quantum  mechanically  so  that  the  limitations  on  the  ultimate  po- 1 
sition  resolution  due  to  atomic  motion  and  wave-mechanical  diffraction  are  included  in  the  analysis. 
The  techniques  utilize  a  miniaturized  form  of  Raman-induced  resonance  imaging  in  which  optical 
fields  are  used  to  make  transitions  from  a  long-lived  initial  state  to  a  long-lived  final  state.  The  final 
state  is  shifted  by  the  spatially  varying  potential  of  an  applied  force  field  in  order  to  correlate  the 
atomic  position  with  its  resonance  frequency.  Spatially  varying  level  shifts  are  obtainable  by  using 
very  large  Zeeman  field  gradients  or  spatially  varying  light  shifts  in  a  small  interaction  volume. 
This  results  in  very  high  spatial  resolution.  The  atomic  transit  time  across  the  optical-field  region  is 
limited  by  focusing  to  an  ideal  diameter  that  minimizes  the  spatial  resolution  length.  The  results  of 
the  analysis  show  that  nanometer  spatial  resolution  of  the  initial-state  position  distribution  is  attain¬ 
able.  Under  appropriate  conditions,  the  final-state  spatial  wave  function  can  take  the  form  of  a 
minimum-uncertainty  Gaussian  wave  packet  obeying  AxAp  =fi/2.  Such  states  may  prove  useful  in 
studying  one-dimensional  wave-packet  motion  in  applied  potentials. 


1.  INTRODUCTION 

Recently  there  has  been  substantial  progress  and  in¬ 
terest  in  the  creation  of  atomic  distributions  which  vary 
over  small  length  scales.  In  atomic  beams,  such  distribu¬ 
tions  arise  by  diffraction,*  periodic  spatial  modulation,^ 
channeling,^  focusing,^  and  cooling.’  In  general,  the  dis¬ 
tributions  exhibit  momentum-space  coherence  and  are  of 
both  practical  and  fundamental  interest.  These  develop¬ 
ments  have  been  spurred  on,  in  part,  by  the  advent  of 
practical  laser-cooling  techniques.^  Important  applica¬ 
tions  include  atomic  interferometry,  gyroscopes,  and  the 
creation  of  submicrometer  structures  by  atomic  deposi¬ 
tion. 

While  techniques  exist  to  create  suboptical-wavelength 
atomic  spatial  modulation  and  interference  in  atomic 
beams,  methods  for  detection  have  been  limited  princi¬ 
pally  to  hot  wires.  Only  in  Ref.  3  is  optical  absorption 
used  to  determine  the  atomic  position  distribution.  In 
that  work,  atoms  are  channeled  in  an  off-resonant  optical 
standing  wave.  Due  to  the  spatial  variation  of  the  light 
.shifts  in  the  standing  wave,  atoms  at  a  node  resonate  at  a 
different  frequency  than  atoms  at  a  peak  in  the  optical 
field.  Thus,  a  two-peak  absorption  spectrum  is  obtained 
in  which  the  height  and  shape  of  the  peaks  vary  when  the 
atoms  are  redistributed  by  channeling,  compared  to  when 
they  are  not  channeled.  Since  the  channeling  field  and 
the  spatially  varying  light  shift  which  causes  the  channel¬ 
ing  vary  with  the  same  scale  of  length,  the  true  periodic 
structure  of  the  channeled  atomic  distribution  could  not 
be  directly  measured  in  the  experiments.  While  the  posi¬ 
tion  resolution  obtained  was  adequate  (although  not 
specified  in  Ref.  3),  it  was  limited  by  the  spontaneous 
linewidth  for  the  transition  employed  as  well  as  by 
Doppler  broadening.  Further,  the  diameter  of  the  in¬ 
teraction  region  (4.6  mm)  was  chosen  to  achieve  channel¬ 


ing  without  heating.  In  this  case,  the  interaction  time  is 
long  enough  for  the  atoms  to  be  redistributed  over  dis¬ 
tances  comparable  to  an  optical  wavelength.  For  op¬ 
timum  position  resolution,  our  calculations  show  that  the 
interaction  time  of  the  atoms  with  the  measurement  field 
must  be  much  shorter,  since  the  atoms  must  move  no 
more  than  the  optimum  resolution  length.  This  is  much 
shorter  than  an  optical  wavelength  for  the  large  light 
shifts  employed  in  the  experiments. 

The  purpose  of  this  paper  is  to  analyze  a  recently  sug¬ 
gested  class  of  techniques’  which  may  prove  useful  in 
achieving  much  higher  spatial  resolution  of  moving 
atoms  using  optical  fields  than  has  been  obtained  previ¬ 
ously.  We  begin  by  reviewing  the  principal  features  of 
the  method  using  heuristic  arguments  and  then  proceed 
to  the  detailed  derivation. 

The  basic  method,  Fig.  1,  is  to  use  a  spatially  varying 
potential  K(x),  which  makes  the  resonance  frequency  of 
simple  atoms  position  dependent,  analagous  to  magnetic 
resonance  (MR)  imaging.  Unlike  conventional  MR  imag¬ 
ing,  the  proposed  method  employs  Raman-induced  reso¬ 
nance  imaging  using  optical  fields.  The  use  of  optical 
fields  permits  the  measurements  to  be  made  in  very  small 
volumes  with  very  large  spatially  varying  potentials, 
which  leads  to  quantum-mechanical  limits  on  the  ulti¬ 
mate  measurement  accuracy.  As  shown  in  Fig.  1,  the 
atoms  enter  the  Raman  interaction  region  in  a  long-lived 
initial  state  (i)  which  is  not  shifted  by  the  potential  Vix) 
and  make  a  transition  via  two  optical  fields  to  a  long- 
lived  final  state  (/)  which  is  shifted  by  the  applied  poten¬ 
tial.  In  this  way,  the  atomic  position  and  its  resonance 
frequency  are  correlated.  The  total  number  of  atoms  em¬ 
erging  in  the  final  state  is  measured  downstream  from  the 
Raman  region.  This  signal  is  proportional  to  the  number 
of  atoms  in  the  initial  state  which  make  a  transition  to 
the  final  state  near  the  position  where  the  optical  fields 
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QUANTUM  THEORY  OF  ATOMIC  POSITION  MEASUREMENT  . . . 
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FIG.  1.  Atomic-position  measurement  by  Raman-induced 
resonance  imaging.  Atoms  enter  the  Raman  interaction  region 
in  an  initial  state  i  with  some  spatial  distribution  along  the  mea¬ 
surement  axis  X  which  is  to  be  measured.  The  initial-state  am¬ 
plitude  Aj(y,x)  is  presumed  to  be  specified  at  the  plane  y  =>>1 
centered  in  the  Raman  fields  fl,,!!;.  A  spatially  varying  poten¬ 
tial  V{x)  shifts  the  final  state  /  without  affecting  the  initial  state 
i  and  correlates  the  Raman  resonance  frequency  with  the  atom¬ 
ic  position  along  the  x  axis.  Atoms  which  are  detected  down¬ 
stream  in  state  /  are  those  which  have  made  the  /  — * /  transition 
in  a  narrow  spatial  region  Ax  where  the  Raman  fields  are  reso¬ 
nant  for  the  I  — */  transition. 


are  resonant  with  the  Raman  transition.  The  goal  is  to 
measure  prepared  incident  atomic  distributions  without 
modification  due  to  trapping  or  channeling,  by  appropri¬ 
ately  limiting  the  interaction  time.  These  techniques  ex¬ 
amine  interesting  features  of  quantum-mechanical  posi¬ 
tion  measurement  and  should  make  it  possible  to  localize 
atoms  to  much  better  than  an  optical  wavelength  for  ap¬ 
plication  to  a  variety  of  interesting  experiments.  For  the 

’^Sm  or  ’’^Yb  systems  we  are  currently  studying,  it 
should  be  possible  to  achieve  position  resolution  of  ==70 
A. 

An  important  feature  of  the  analysis  is  that  it  deter¬ 
mines  the  ultimate  spatial  resolution  attainable  for  mov¬ 
ing  atoms  in  large  spatially  varying  level  shift  gradients 
which  correlate  the  atomic  resonance  frequency  with  po¬ 
sition.  The  results  should  have  important  bearing  on  ex¬ 
periments  which  measure  the  interaction  of  atoms  with 
surfaces.* 

Finally,  it  is  shown  below  that  the  center-of-mass  wave 
function  for  atoms  which  emerge  in  the  final  state  is  high¬ 
ly  localized  and,  under  appropriate  conditions,  can  take 
the  form  of  a  minimum  uncertainty  Gaussian  wave  pack¬ 
et.  These  results  complement  recent  work  on  atomic  and 
molecular  internal  state  wave  packets.'’ 

For  optimum  position  resolution,  it  is  necessary  that 
both  radiative  and  Doppler  broadening,  as  well  as  the 
laser  linewidth,  be  negligible  compared  to  transit-time 
broadening.  In  some  cases,  a  one-photon  transition  from 
the  ground  state  is  adequate.  This  would  require  a  long- 
lived  excited  State,  large  spatially  varying  level  shifts,  and 
high  atomic-beam  collimation.  Alternatively,  a  Raman- 
induced  resonance  between  long-lived  ground  states  can 
be  used,’  as  mentioned  above  and  described  below. 

To  explain  the  Raman-induced  resonance  imaging 
method  in  more  detail,  suppose  that  an  atomic  beam  is 
prepared  in  an  initial  long-lived  internal  state  (i)  with 


some  interesting  spatial  distribution  perpendicular  to  the 
atomic-beam  propagation  axis.  To  measure  the  spatial 
distribution,  the  atoms  enter  an  interaction  region  con¬ 
taining  two  optical  fields  of  frequencies  ni,fl2  which  can 
induce  Raman  transitions  between  the  long-lived  initial 
state  and  a  long-lived  final  state  /.  The  interaction  region 
also  contains  a  spatially  varying  potential  which  makes 
the  Raman  resonance  frequency  dependent  on  the  spatial 
position  of  the  atom.  For  simplicity,  suppose  that  the 
spatially  varying  potential  is  due  to  a  uniform  force  F 
which  affects  only  the  final  state  ( /)  of  the  atom.  The  po¬ 
tential  is  taken  to  vary  linearly  along  the  x  axis,  which  is 
perpendicular  to  the  atomic-beam  propagation  direction, 
and  which  ultimately  will  become  the  measurement  axis. 
For  concreteness,  let  XpFbe  the  shift  of  the  final  state  due 
to  a  uniform  Zeeman  field,  and  —  xF  be  the  shift  of  the 
final  state  due  to  an  additional  Zeeman  shift  or  light  shift 
which  varies  linearly  in  space.  The  net  potential  can  be 
written  as 

F(x  — Xo)=  — (x  — XqIF  .  (1) 

Further,  suppose  that  the  atoms  which  enter  this  poten¬ 
tial  in  the  internal  state  (i)  are  not  shifted.  In  this  case, 
the  i-f  resonance  frequency  varies  linearly  in  space,  so 
that  a  resonance  field  of  sufficient  frequency  resolution 
can  cause  i  — *•/  transitions  in  a  highly  localized  region  Ax 
along  the  x  axis.  The  position  resolution  Ax  is  deter¬ 
mined  by  the  potential-energy  change  AK=FAx  and  the 
frequency  resolution  Aw  according  to 

•^Ax=Aw=~  ,  (2) 

n  T 

where  T*  is  the  effective  interaction  time  which  deter¬ 
mines  the  frequency  re.solution.  If  a  one-photon  optically 
allowed  transition  is  employed,  then  the  frequency  resolu¬ 
tion  will  be  limited  by  Doppler  and  radiative  broadening 
as  well  as  by  laser-frequency  jitter.  All  of  these 
deficiencies  are  remedied  by  using  a  copropagating-wave 
Raman-induced  resonance  between  long-lived  states  (i) 
and  (/)  via  some  off-resonant  intermediate  state  (/)  as 
shown  in  Fig.  2. 

The  optical  fields  fli.flj  are  taken  to  copropagate.  By 
choice  of  the  difference  frequency  flj  — fl^,  a  Raman  reso¬ 
nance  with  the  transition  frequency  can  be 

made  to  occur  at  the  point  x  =Xo  along  the  measurement 
axis.  The  Raman  fields  are  derived  by  A  /O  modulation 
from  the  laser  field  <u  ,so  that  the  difference  frequency  is 
easily  made  very  stable.  By  u.sing  an  off-resonant  Raman 
transition  and  copropagating  fields,  the  laser  jitter,  beam 
divergence,  Doppler  shifts,  etc.,  drop  out,  since  only  the 
difference  frequency  enters.  This  technique  was  em¬ 
ployed  previously  to  obtain  Raman-induced  Ramsey  res¬ 
onances  using  optical  fields,  where  100-Hz  resolution  was 
readily  obtained.'”  By  monitoring  the  total  count  rate  of 
atoms  in  state  /  emerging  from  the  Raman  region  versus 
uniform  Zeeman  shift  (which  varies  xj,),  the  incident  spa¬ 
tial  distribution  in  state  /  is  determined.  Only  initial  state 
atoms  near  x  =Xo,  the  point  at  which  n,  — n2  =  tu^,,  con¬ 
tribute  to  the  signal. 

The  use  of  optical  fields  in  the  measurement  leads  to  a 
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FIG.  2.  Level  scheme  and  Raman  fields. 


number  of  new  features  which  are  li.stcd  briefly  here. 

(1)  Spatially  Varying  Level  Shifts.  Two-wire  magnetic- 
field  gradients  as  developed  for  atom-deflection  applica¬ 
tions  can  be  used  to  achieve  spatially  varying  level  shifts 
F/h  — 10'°  Hz/cm  over  millimeter  length  scales.  Over 
smaller  maximum  distances,  very  large  spatially  varying 
level  shifts  can  be  obtained  by  using  the  light  shifts  of  in¬ 
tense,  nonresonant  intersecting  laser  beams  {not  the  Ra¬ 
man  beams).  In  this  case,  more  than  lO'^  Hz/cm  is  feasi¬ 
ble.  Spatially  varying  shifts  this  large  have  been  used  in 
atom-channeling  experiments,  but  application  to  position 
measurements  which  yield  the  incident  spatial  distribu¬ 
tion  have  not  been  explored.  When  the  interaction  time 
is  optimally  limited  as  described  below,  exceedingly  high 
spatial  resolution  is  attainable. 

(2)  Spatial  Resolution.  When  long-lived  states  (/,/)  are 
chosen  for  measurement,  the  interaction  time  T*  is  deter¬ 
mined  by  the  atom  transit  time  T  across  the  Raman  re¬ 
gion.  In  this  case,  the  use  of  optical  fields  allows  the  in¬ 
teraction  time  to  be  limited  to  an  optimum  value  by 
focusing.  The  optimum  value  of  T  is  determined  by  the 
criterion  that  the  atom  not  move  across  the  ideal  resolu¬ 
tion  length  during  the  measurement.  The  simplest  limit 
is  due  to  the  atomic  velocity  component  along  the 
measurement  axis.  In  this  case,  the  optimum  interaction 
time  is  T  =  Ax/v^.  With  Eq.  (2),  the  optimurn  spatial 
resolution  is  given  by 

Ax„|  =  /F  .  (3) 

For  modest  atomic-beam  collimation  and  a  two-wire 
magnetic-field  gradient,  suboptical-wavelength  position 
resolution  is  readily  attainable. 

For  a  highly  collimated  atomic  beam,  Uj,  no  longer  lim¬ 
its  the  interaction  time.  In  this  case,  the  acceleration  due 
to  the  applied  force  and  diffraction  out  of  the  resolution 
length  limit  the  smallest  regions  which  are  measurable. 
For  the  acceleration  limit,  the  optimum  value  of  the  in¬ 
teraction  time  is  such  that  {F/2M)T^  =  ^x.  For  the 
diffraction  limit,  the  interaction  time  must  be  such  that  a 
particle  localized  into  a  region  of  length  Ax  spread  only 
over  distances  comparable  to  Ax,  i.e.,  (h/lM tsx)T  —  S.x. 
In  either  case,  with  Eq.  (2),  the  optimum  resolution  is 


which  we  will  refer  to  as  the  acceleration  limit.  To 
achieve  this  limit,  we  must  have  v^  <ff/(2Aif  Axj„ci  ** 
shown  below.  For  a  lO'^-Hz/cm  light-shift  gradient,  the 


optimum  resolution  is  about  70  A  for  samarium  and  the 
required  collimation  for  a  beam  of  speed  6X  10^  cm/s  is 
SO  /irad,  which  is  not  too  stringent. 

(3)  Uncertainly  limit.  Since  the  force  due  to  the  spatial¬ 
ly  varying  potential  acts  on  the  atom  only  after  it  makes  a 
transition  from  the  initial  unshifted  state,  the  momentum 
uncertainty  imparted  by  the  force,  Ap,  is  approximately 
equal  to  the  maximum  momentum  imparted,  FT.  Using 
Eq.  (2)  we  see  that  the  final  state  is  created  with  a  spatial 
resolution  and  momentum  spread  which  obey  the  uncer¬ 
tainty  relation  AxAp^ff.  In  Sec.  IIIB  it  is  shown  that 
minimum  uncertainty  packets  can  be  created  under  ap¬ 
propriate  conditions. 

(4)  Detection.  The  method  for  detecting  the  final  state 
( /)  need  not  have  any  spatial  resolution  since  the  mea¬ 
surement  is  done  by  the  Raman  region.  Only  the  spatial 
integral  of  the  final-state  flux  is  needed,  therefore  a 
variety  of  final-state  detection  schemes  are  possible. 
These  include  resonance  fluorescence,  and  single-atom 
methods  such  as  photon  burst  and  photoionization, 
which  have  been  successfully  applied  to  trace-element 
measurement. 

II.  THEORY 

In  the  following,  we  consider  the  Raman-induced  reso¬ 
nance  imaging  scheme’  in  detail.  Schemes  based  on  one- 
photon  transitions,  when  appropriate,  yield  analagous  re¬ 
sults  and  need  not  be  derived  in  depth.  A  three-level 
atom  with  states  /,/,/,  as  shown  in  Fig.  2,  enters  a  region 
with  a  spatially  varying  potential  which  exerts  a  force  on 
the  state  /  only.  The  atom  enters  the  region  in  state  i 
which  does  not  interact  with  the  spatially  varying  poten¬ 
tial  and  makes  a  Raman  transition  via  two  off-resonant 
light  fields  to  the  final  state.  The  transition  occurs  in  a 
localized  region  determined  by  the  resonance  condition 
that  the  initial-  lo-final-state  frequency  oiy,  equal  the 
difference  frequency  of  the  two  applied  optica)  fields.  The 
total  (spatially  integrated)  final-state  flux  emerging  from 
the  Raman  region  is  monitored  downstream  as  a  function 
of  the  Raman-field  difference  frequency  or  as  a  function 
of  (Off  which  can  be  tuned  by  means  of  an  applied  uni¬ 
form  magnetic  field. 

In  this  section,  the  integrated  flux  of  atoms  in  state  f 
exiting  the  Raman  region  is  calculated  at  distances  large 
compared  to  the  diameter  of  the  Raman  regions  which 
are  assumed  to  be  Gaussian.  The  steady-state  amplitude 
for  the  final  state  is  calculated  first,  where  it  is  assumed 
that  the  wave-mechanical  motion  is  important  for  the  x 
motion.  An  eikonal  approximation  is  assumed  to  be  ade¬ 
quate  for  the  y  coordinate  along  the  atomic-beam  axis. 
The  Raman  transition  frequency  is  taken  to  be  very 
small,  so  that  the  difference  in  the  wave  vectors  of  the  op¬ 
tical  fields  can  be  taken  as  zero.  In  this  case,  the  net 
recoil  momentum  of  the  atom  in  making  the  Raman  tran¬ 
sition  is  negligible.  Finally,  the  optical  fields  are  assumed 
to  be  far  enough  off  resonance  with  the  excited  state  that 
the  excited-state  amplitude  can  be  adiabatically  eliminat¬ 
ed  from  the  coupled  amplitude  equations  which  are  de¬ 
rived  below. 

The  Hamiltonian  for  the  three-level  system  of  Fig.  2  is 
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taken  as 


2jVi 


where  //**”  is  the  internal  Hamiltonian,  P  is  the  atomic 
center-of-mass  momentum  operator,  is  the  external 
force  operator,  and  x  is  the  position  along  the  measure¬ 
ment  axis.  The  electric  dipole  operator  is  fi  and  the  ap¬ 
plied  field  is  of  the  form 


E(R,f)=  +C.C. 


+  --e 


-1-c.c.  e 


-\/2I(y-y,)/dH‘ 


(6) 


where  the  Raman  fields  are  assumed  to  copropagate  and 
the  diiference  frequency  n,  — 112  is  assumed  sufficiently 
small  that  the  wave  vectors  are  taken  to  be  of  equal  mag¬ 
nitude.  The  polarizations  Cj  2  are  chosen  so  that  each 
field  interacts  with  only  one  transition.  The  Raman-field 
region  is  taken  to  be  centered  at  a  position  y,  along  the 
atomic-beam  axis  and  the  laser  1  /e  intensity  radius  along 
the  y  axis  is  d.  Variation  in  the  fields  along  directions 
perpendicular  to  the  atomic-beam  axis  is  assumed  negligi¬ 
ble.  The  atomic  wave  function  is  taken  to  be  of  the  form 

n  ^  ^ 

a) 


where  the  velocity  of  the  particular  atom  being  con¬ 
sidered  is  v—Hk/M  prior  to  interaction  and  it  is  assumed 
that  an  initial  velocity  average  will  be  performed  at  the 
end  of  the  calculation.  The  internal  state  of  energy 
is  |n>.  With  these  assumptions,  the  time-dependent 
Schrodinger  equation  yields  equations  of  motion  for  the 
state  amplitudes  A„  given  by 

m 

where  the  force  F„  =Fb„f  acts  only  on  the  final  state. 
We  take  the  Raman  difference  frequency  equal  to  the 
/■  — ►/  frequency  in  the  absence  of  the  force:  fli  — n2=<Uy/. 
Then  we  seek  solutions  for  which  the  initial-  and  final- 
state  amplitudes  are  time  independent  but  vary  in  space. 
It  is  assumed  that  the  velocity  perpendicular  to  the  laser 
beams  is  large  (i.e.,  a  thermal  speed),  while  the  velocity 
u,  along  the  measurement  axis  x  is  small.  In  this  case,  it 
is  convenient  to  make  an  eikonal  approximation  for  they 
motion,  retaining  only  the  first  spatial  derivative  term, 
while  for  the  x  direction,  both  the  first  and  second  x 
derivatives  are  retained  to  take  into  account  packet 
.spreading.  Equation  (8)  is  then  easily  solved  in  perturba¬ 
tion  theory  by  finding  the  Green’s  function  for  states  i 


and  /  and  treating  the  right-hand  side  as  a  source. 

To  proceed,  Eq.  (8)  is  rewritten  explicitly  using  the 
rotating-wave  approximation  and  takes  the  form 

fl'ifR)  -/q.R-i(„  -n,)( 

= - "  '  /l,(R,f), 


_  —  iq-R  — 


'  A,(R,t}  ,  (9) 


3  .  _  ifi 

— -fvV-— — — r  R,t) 

dt  2M  dx^ 


2i 


•^,(R.f) , 


where  ^l|i{R)=^l,j^t^S^{R)/fi  and  fl/^-lR) 
=tt,f-C2Si{R)/fi  and  the  spatial  dependence  of  the  fields 
is  given  by  the  Gaussian  distribution  of  Eq.  (6).  The  con¬ 
vective  derivative  term  describes  a  straight-line  trajectory 
for  constant  velocity,  the  second  derivative  term  in  x  al¬ 
lows  for  diffraction  (wave-packet  spreading),  and  the 
force  term  takes  into  account  the  spatially  varying  poten¬ 
tial  which  affects  only  the  final  state. 

It  is  assumed  that  the  Raman  fields  are  detuned  far 
from  resonance  with  the  intermediate  state  so  that 
coherent  transitions  from  the  initial  state  to  the  final  state 
arc  dominant  over  incoherent  optical  pumping.  In  this 
case,  the  intermediate-state  amplitude  A  /  can  be  adiabat- 
ically  eliminated  from  Eqs.  (9).  Suppose  that  fl,- 
f22  that  o) fli — ^ If  fl2  =  A.  In 

the  large-detuning  limit,  we  have  A»y,,  q-v,  fl,,, 
etc.  Letting  Af{t)=a,{t)txp{~i^t)  in  Eqs.  (9),  the 
intermediate-state  amplitude  is  approximately  given  by 


With  Eq.  (10),  the  evolution  equations  (9)  for  the  initial- 
and  final-state  amplitudes  take  the  form 
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Equations  (11)  are  readily  solved  in  perturbation 


theory.  This  yields  a  spatial  resolution  function  (see 
below)  which  includes  the  effects  of  both  classical  motion 
and  diffraction  of  the  moving  atoms.  In  lowest  order,  the 
initial-state  amplitude  Aj^Hy.x}  can  be  taken  to  propa¬ 
gate  neglecting  the  Raman  fields.  Further,  the  velocity 
will  be  taken  to  have  only  x  and  y  components.  It  is  as¬ 
sumed  that  the  atomic  distribution  does  not  vary  in  the  z 
direction.  Then,  the  initial-state  amplitude  satisfies 

A.  J _ i±JL 

"’^dy  ^^dx  2M  dx^ 

The  solution  to  this  equation  is  determined  by  free- 
particle  propagation,  where  y/v^  plays  the  role  of  the 
time.  If  the  amplitude  A,  is  specified  in  the  plane  located 
at  y ,  along  the  atomic-beam  axis,  then  the  amplitude  at 
an  arbitrary  position >>2. .*2  evidently  given  by 


Aj°\y,x)=Q  .  (12) 


J 


f^^dx\Ai{y^,x\)  f  dkf 


I  ~-'l  “"1*^2  -yj  -/(A*  J/2/nl(y2  -y, 

2ir 


(13) 


To  lowest  order  in  the  Raman  fields,  the  term  on  the  right-hand  side  of  the  A^  evolution  equation  (11)  proportional 
to  A^  may  be  neglected,  since  A^  is  assumed  to  be  initially  zero,  prior  to  a  Raman  transition  from  state  /.  (That  is,  we 
neglect  the  light  shift  and  optical  pumping  due  to  the  weak  nonresonant  Raman  fields.  In  practice,  the  Raman-field  in¬ 
tensities  are  adjusted  to  give  equal  light  shifts  which  cancel  for  the  1  — ►/  transition.)  The  solution  for  the  final-state  am¬ 
plitude  in  the3»,x  plane  is  then  readily  obtained  by  means  of  the  Green’s  function  which  satisfies 
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The  final-state  amplitude  is  given  by 


4 

Ys  .. 

- lA 

2 

(15) 


where  the  Rabi  frequencies  are  taken  to  be  Gaussian  distributions  centered  at  >>,  along  the  y  axis  [see  Eq.  (6)]  but  in¬ 
dependent  of  X  for  simplicity  (i.e.,  the  measurement  length  is  taken  to  be  small  compared  to  the  radius  and  the  Rayleigh 
length  of  the  Raman  beams). 

As  shown  in  Appendix  A,  the  required  Green’s  function  for  propagating  in  a  uniform  force  field  is  obtainable  by 
transforming  to  momentum  space  and  is  given  by 


Gf{y,y2;x,x2)= 


eiy  -y2 ) 
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dk~  exp 
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(16) 


In  this  form,  the  Green’s  function  has  a  simple  physical  interpretation.  Neglecting  wave-packet  spreading,  integration 
over  k  for  the  first  exponential  factor  yields  a  5  function  requiring  the  atom  to  follow  the  classical  trajectory  during  the 
time  {y  — >’2  )/Uy.  The  second  exponential  takes  wave-packet  spreading  into  account.  The  last  exponential  factor  is  just 
the  phase  acquired  by  the  final  state  of  the  atom  in  propagating  through  the  classical  path  in  the  potential  due  to  the 
force  field. 

In  Appendix  B,  the  integrated  final-state  flux  emerging  from  the  Raman  region  is  evaluated  starting  from  Eq.  (15)  and 
yields  the  following  result; 
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(17) 


Equation  (17)  gives  the  general  result  for  the  integrated 
flux  in  the  final  state  in  terms  of  an  arbitrary  initial-state 
amplitude  specified  at  the  center  of  the  Raman  region 
y  =>'|.  Ajiy^fX ).  To  obtain  the  final-state  count  rate,  this 
result  must  be  multiplied  by  the  atomic  density,  the 
atomic-beam  velocity,  and  the  transverse  dimension  of 
the  atomic  beam  perpendicular  to  the  measurement  axis. 

III.  RESULTS 

In  the  following  sections,  the  spatial  resolution  ob¬ 
tained  with  the  Raman-induced  resonance  imaging 
method  is  examined  first,  in  order  to  verify  the  statements 
made  in  the  Introduction  based  on  heuristic  arguments. 
Then,  the  final-state  wave  function  is  examined  to  deter¬ 
mine  the  momentum  and  position  distribution  of  the 
final-state  atoms.  For  a  plane-wave  input  in  the  initial 
state,  it  is  shown  that  under  appropriate  circumstances,  a 
minimum-uncertainty  final-state  wave  packet  may  be  ob¬ 
tained. 


A.  The  spatial  resolution  function 

Using  the  result  given  by  Eq.  (17),  the  spatially  In¬ 
tegrated  final-state  flux  can  be  determined  for  any  given 
spatially  varying  initial-state  amplitude  Ai{yi,x].  This 
determines  the  dependence  of  the  signal  on  the  applied 

- _ ) 
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uniform  magnetic  field  (i.e.,  on  the  resonant  point  Xq). 
Further,  in  an  important  special  case,  described  below, 
this  result  is  to  determine  the  general  form  of  the  sam¬ 
pling  function  which  measures  the  incident  atomic  distri¬ 
bution  in  state  We  will  refer  to  this  sampling  function 
as  the  resolution  function. 

The  Gaussian  factors  in  the  integrand  of  Eq.  (17)  con¬ 
strain  u  ^  to  distances  of  order  d.  If  the  initial-state  am¬ 
plitude  does  not  vary  over  distances  of  order 
{F/2M){d /Vy)^,  then  the  dependence  of  the  amplitudes 
on  Uj.  may  be  neglected.  This  requires  that  the  interac¬ 
tion  time  d/Vy  in  the  Raman  fields  be  not  too  large,  so 
that  the  final  state  does  not  accelerate  over  distances 
comparable  to  the  length  scale  over  which  the  initial- 
state  amplitude  varies.  In  this  case,  the  u  +  integration  is 
readily  done  and  the  integrated  flux  then  takes  the  form 

dx\Af{y,x)\^  = 

X  f  dx\R\[x  \  — Xq) 

xM,(y,,x;)(^  (18) 

where  the  resolution  function  R ,  is  given  by 
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Physically,  the  meaning  of  this  result  is  clear.  The  ex¬ 
ponent  in  the  Gaussian  factor  decreases  as  the  Raman 
beam  radius  d  is  made  larger  and  the  interaction  time 
longer,  corresponding  to  higher-frequency  resolution. 
This  tends  to  increase  the  spatial  resolution  and  narrow 
the  resolution  function.  However,  as  d  is  increased,  the 
motion  along  the  measurement  axis  due  to  tends  to  in¬ 
crease  the  exponent  and  broadens  the  resolution  function. 
The  acceleration  of  the  final  state  leads  to  a  phase  factor 
in  the  integrand  which  is  cubic  in  the  time  (u_  /Vy),  as 
one  would  expect  using  classical  arguments.  In  addition. 


diffraction  also  leads  to  a  cubic  time  dependence  of  the 
phase. 

This  result  can  be  cast  in  a  simpler  form  by  using  the 
natural  length  scales  for  this  problem  which  are  discussed 
in  the  Introduction.  The  form  of  the  resulting  resolution 
function  serves  to  justify  the  heuristic  treatment  of  the 
resolution  limits. 

It  is  convenient  to  express  distance  along  the  measure¬ 
ment  axis  in  units  of  the  acceleration  limited-resolu¬ 
tion  length,  Eq.  (4).  The  resolution  length  in  turn  leads 
to  the  natural  unit  of  transit  time  across  the  Raman 
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region,  Tq.  Using  {F/2M)TI  =  Ax^^^^^  one  obtains 
7’o  =  (2A//i/F^)'^\  The  corresponding  unit  of  distance 
along  the  atomic-beam  axis,  is  just  Defining 

u_  —r]do,  and  expressing  x\  — Xq  in  units  of  the 

resolution  function  can  be  recast  into  the  following  form: 

1  x'l  — X(i 

/IqI-*  1 ^  XT  > 


j{s)=  f  ”  t/n— -isgn(ns>)-i(T(V6) 

-  00  2n 

=  r  dn— e  “"'’'^'’^'cos  sgn{F)sr]+  ^ -  ,  (21) 

•'  n  IT  6 

with  sgn(F)=  +  1(  -  1 )  for  a  positive  (negative)  force. 
Note  that  J dx  |/?n(.x',  — Xn)=  1.  The  parameters  appear¬ 
ing  in  Eq.  (21)  are  given  by 

do  d'  4 
a  =  — rH - rr  , 

d^  dl 


Ax„,,= 


do  =  Vy 


In  terms  of  Fq,  the  integrated  final-state  flux  takes  the 
form 


J"  dxlA^ijr.x) 
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-4,(^1, x'=Xo)|^  ,  (23) 


where  the  approximate  result  is  valid  if  the  resolut'on 
function  is  narrow  compared  to  the  length  scales  over 
which  the  atomic  distribution  varies.  The  Raman  transi¬ 
tion  probability  is  given  by 


Note  that  the  factor  which  multiplies  the  Raman  transi¬ 
tion  probability  in  Eq.  (23)  is  of  the  order  of  the  position 
resolution.  The  approximate  result  could,  of  course,  have 
been  written  down  on  the  basis  of  physical  arguments, 
without  detailed  calculation.  The  important  feature  is 


that  the  detailed  shape  of  the  resolution  function  is  ob¬ 
tained. 

According  to  Eqs.  (20)  and  (21),  the  resolution  function 
will  be  approximately  Gaussian  if  the  param'^ter  a  is 
large  compared  to  unity,  so  that  the  cubic  dependence 
of  the  phase  on  t}  can  be  neglected.  When  the 
velocity  is  so  large  that  Ax„|  >  Ax  which  implies 

></2MAx,^„i,  the  factor  r*  in  a  of  Eq.  (22)  quickly 
becomes  large  and  Gaussian  dependence  dominates  the 
integrand.  In  this  limit,  the  resolution  function  takes  the 
form 


/(s)  = 


'Z  27ra 


(25) 


According  to  Eq.  (22),  the  parameter  a  is  minimized  for  a 
certain  value  of  the  Raman-region  diameter  d.  Physical¬ 
ly,  the  frequency  resolution  increases  as  d  is  increased. 
This  tends  to  increase  the  position  resolution.  However, 
the  velocity  along  the  measurement  axis  causes  the  po¬ 
sition  to  move  during  the  measurement.  Hence,  an  op¬ 
timum  diameter  exists.  Differentiating  a  with  respect  to 
dy  and  setting  the  result  equal  to  zero,  one  obtains 
dap\~d(i/r  and  a„p,  =  2rT  The  corresponding  resolution 
function  then  is  given  by  Eqs.  (20)  and  (25)  as 


Rq{x  ^  Xq  ) 


Ax^i^tt 


(26) 


and  d„p,=(u^/Uj,)Ax„i-  Physically,  the  atom  travels  a 
distance  along  the  measurement  axis  just  equal  to  the  op¬ 
timum  resolution  length  during  the  transit  time  across 
the  Raman  region.  Note  that  the  neglect  of 
Fu  +  M  _  /2MVy  in  the  argument  of  the  initial  amplitude 
Aj  in  Eq.  (17)  is  justified  if  A,  varies  slowly  compared  to 
Fdlpt/2MVy=fi/2Mi,  <  Ax,„,|  in  this  limit. 

In  the  opposite  limit  when  Ax„f|  <  Ax,^„i,  i.e.,  when 
Uj,  <^/2Af  Ax,„j|,  the  ratio  r  is  less  than  unity  and  the 
parameter  a^d^/d^.  Hence,  as  the  Raman-region  di¬ 
ameter  d  is  increased  to  order  d^  or  greater,  a  decreases, 
and  the  functions  Us)  and  narrow.  As  a— *0,  it  is 
clear  from  the  exact  equation  (21)  for  /(s)  that  the  natu¬ 
ral  length  scale  for  s  is  unity,  since  the  cubic  phase  plays 
the  dominant  role  in  the  integrand.  This  corresponds  to 
spatial  resolution  comparable  to  the  acceleration  length. 
The  function  I{s)  is  plotted  in  Fip  3  for  various  values  of 
a. 

When  a  is  of  order  unity  (1-2),  the  shape  is  nicely  lo¬ 
calized  and  is  nearly  Gaussian.  However,  if  a  is  made 
small,  corresponding  to  a  large  Raman-region  diameter, 
the  function  /(j)  is  asymmetric.  It  oscillates  on  one  side 
and  decays  on  the  other.  Physically,  this  is  linked  to  the 
fact  that  for  a  positive  force,  the  final-state  amplitude 
created  at  x'  to  the  left  of  the  resonance  position  Xq 
within  a  distance  of  order  Ax^^ei  accelerates  to  the  right 
and  interferes  with  the  amplitude  created  at  Xg.  Hence, 
the  resolution  function  oscillates  on  the  left  side  of  Xq  for 
a  positive  force.  However,  final-state  amplitude  created 
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for  x’  to  the  right  of  Xq  accelerates  further  to  the  right 
away  from  jcq  and  hence  does  not  interfere.  Thus,  the 
resolution  function  rapidly  dies  off  on  this  side  of  Xg  for  a 
positive  force.  This  behavior  is,  of  course,  related  to  the 
wave  function  for  penetration  of  a  potential  barrier  which 
varies  linearly  in  position,  but  is  opposite  to  what  one 
might  expect  because  the  final  state  is  accelerated  rather 
than  the  initial  state. 

B.  Minimum-uncertainty  spatiai  wave  packets 

The  process  of  measuring  the  initial-state  distribution 
via  the  Raman-induced  resonance  imaging  method  natu¬ 
rally  creates  the  final-state  spatial  wave  function  in  the 
form  of  a  localized  packet,  since  the  i~*f  transition 
occurs  in  a  highly  localized  spatial  region.  Hence,  the 
method  is  well  suited  to  the  creation  and  study  of  spatial 
wave  packets  and  it  is  of  interest  to  determine  the 
momentum  and  position  spread  of  such  a  state.  For  this 
purpose,  we  assume  for  simplicity  that  the  initial  state 
takes  the  form  of  a  plane  wave  with  Ai(yy,x}=i.  In  this 
case,  the  final-stats  spatial  wave  function  after  the  Raman 
region  where  y>yi  can  be  evaluated  starting  from  Eq. 
(B6).  The  X integration  is  straightforward  and  yields  the 
unnormalized  state  for  an  atom  of  initial  velocity 
v=^k/Af  as  il)f^{y,x]=e''‘'^A^(y,x).  Choosing  a  nor¬ 
malization  constant  so  that  f  Z^dx\il>^(y,x)\^=\  yields 
the  final-state  wave  function  in  the  form 
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FIG.  3.  Spatial  resolution  function.  The  s  axis  denotes  posi- 
tion  along  the  measurement  axis  in  units  of  the  optimum  ac¬ 
celeration  limited  resolution  length  (see  text). 
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(27) 


Using  Eq.  (27),  the  momentum  spread  of  the  final-state 
wave  packet  is  readily  evaluated,  since  differentiation 
with  respect  to  x  just  differentiates  the  phase  factors  and 
evaluation  of  the  expectation  values  then  requires  only 
Gaussian  integrals.  The  momentum  spread  is  found  to  be 

12 


reasonable  approximation.  With  adequate  beam  collima- 
tion,  the  acceleration  limit  is  attained  and  yields 
(^p)^=[{.Fd)nVy)f. 

The  position  spread  also  can  be  evaluated  from  Eq.  (27) 
using  integration  by  parts  and  the  identity 


(Ap)^=^— p  ^  =  )^  + 


Fd 

2u„ 


(28) 


f  ”  du  u  "  expi 

—  OD 


F  yt-y2 


where  the  first  term  arises  from  a  thermal  average  over 
the  transverse  velocity  spread  and  the  spread  in  the  longi¬ 
tudinal  velocity  Vy  along  the  atomic  beam  is  ignored  for 
simplicity.  For  a  supersonic  atomic  beam,  this  is  a 


2ir 


\F/fiv, 


a 


iF 


yi 


8(^2 -pi ) 


(29) 


This  yields 


fiVy 

~Fd 


+ 


Fd^ 


2Mvl 


P-Pi 


y  -y\ 


Fv;  dHy-y^) 
Mvl  4 


(30) 


where  the  overline  on  the  Uj,  terms  indicates  a  thermal 
average  over  the  initial  transverse  velocity  spi  ead  and  the 
Uy  spread  is  again  ignored.  For  these  equations  to  be  val¬ 
id,  the  observation  point  must  be  after  the  Raman  region. 


so  thatp  — p,  must  be  somewhat  greater  than  d. 

Using  the  position  and  momentum  spread,  the  uncer¬ 
tainty  product  can  be  obtained,  ihe  most  interesting 
case  is  that  of  high  collimation  where  the  spread  in 
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satisfies  Av^  <fi/(M Ax^„i)  and  the  accel' 'ation  limit  is 
achieved.  In  this  case,  the  transverse  coherence  length  of 
the  incident  beam  is  greater  than  and  the  force  F 

imparts  a  momentum  larger  than  the  transverse  thermal 
spread.  Neglecting  the  transverse  velocity  spread  terms, 
the  uncertainty  product  then  takes  the  form 


(Ax)^(Ap)^=-^ 


1  + 


J’-J’l 

2 

d 

(31) 

where  c/q  natural  unit  of  distance  along  the 

atomic-beam  axis  in  the  acceleration  limit,  Eq.  (22).  For 
observation  points  not  too  far  after  the  Raman  region, 
y  — 3>i  is  of  the  order  of  Raman  region  size  d  and  Eq.  (31) 
shows  that  a  minimum  uncertainty  state  is  rapidly  ap¬ 
proached  as  the  radius  d  is  reduced  somewhat  below  d^. 
Physically,  this  is  easy  to  understand.  When  the  radius  d 
of  the  Raman  region  is  reduced  below  the  size  dg  which 
achieves  optimum  position  resolution  in  the  acceleration 
limit,  atomic  acceleration  and  diffraction  during  the  tran¬ 
sit  time  across  the  Raman  region  are  negligible.  In  this 
case,  the  atom  in  its  rest  frame  sees  a  Gaussian  frequency 
pulse  without  chirping.  Due  to  the  linear  spatial  varia¬ 
tion  of  the  final-state  potential,  the  Gaussian  frequency 
distribution  selects  0'’.t  a  Gaussian  spatial  distribution  for 
the  final-state  wave  packet.  This  follows  from  Eq.  (27) 
when  the  quadratic  terms  in  the  phase  due  to  u,  and  cu¬ 
bic  terms  due  to  the  acceleration  are  neglected,  which 
leads  to  a  Gaussian  spatial  wave  function  for  the  final 
state. 


IV,  CONCLUSIONS 

The  calculation  of  the  form  of  the  resolution  function 
and  of  the  final-state  wave  function  above  'justifies  the 
statements  made  in  Ref.  7  on  the  basis  of  heuristic  argu¬ 
ments.  The  results  for  the  ultimate  spatial  resolution 
Ax,^i  show  that  very  high  spatial  resolution  is  attain¬ 
able  for  large  spatially  varying  final-state  level  shifts.  For 
example,  a  samarium  atomic  beam  of  nominal  speed 
6X  10^  cm/s  in  a  level  shift  gradient  of  lO'^  Hz/cm  yields 
a  spatial  resolution  =70  A.  To  obtain  this  resolution,  it 
is  necessary  that  AjCvji  <  Axr„„(  which  is  equivalent  to 
U;,  <fl/(2AfAx,e„|).  For  the  above  conditions,  this  re¬ 
quires  v^<2  cm/s  corresponding  to  a  beam  collimation 
better  than  50  fim.  The  optimum  Raman  beam  intensity 
1/e  radius  d  is  about  140  /xm  in  this  case.  Note  that  the 
Doppler  shifts  due  to  divergence  of  the  Raman  beams 
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cancel  if  the  beam  radii  are  matched  at  the  focus. 

For  large  spatially  varying  level  shift  gradients,  the  fre¬ 
quency  resolution  required  to  obtain  the  optimum  spatial 
resolution  is  not  as  stringent  as  for  small  gradients.  In 
this  case,  a  simple  one-photon  transition  may  be  applic¬ 
able  with  long-lived  excited  states  to  achieve  very  high 
spatial  resolution.  For  example,  in  '^^Yb,  this  method 
can  be  applied  to  excite  the  'Sg— transition  and  ei¬ 
ther  absorption  or  resonance  fluorescence  in  the  interac¬ 
tion  region  can  be  directly  measured  by  modulating  the 
level  shift  gradient  as  done  by  Salomon  et  al.^  In  this 
case,  for  suboptical  wavelength  spatial  resolution  the 
recoil  momentum  for  a  one-photon  process  will  be  negli¬ 
gible  compared  to  the  momentum  imparted  by  the  spa¬ 
tially  varying  potential. 

An  important  application  of  the  Raman-induced  reso¬ 
nance  imaging  method  is  in  the  development  of  two-point 
correlation  techniques.  By  employing  two  pairs  of  Ra¬ 
man  difference  frequencies,  it  is  possible  to  access  two 
separated  spatial  points  with  high  spatial  resolution.  The 
final  states  for  each  of  the  regions  can  be  made  either 
identical  or  different  so  that  the  contribution  of  each  re¬ 
gion  to  the  total  detected  signal  can  be  made  either  dis¬ 
tinguishable  or  indistinguishable.  The  use  of  two-point 
correlation  measurements  opens  up  the  possibility  of 
achieving  the  atomic  analogue  of  the  optical  intensity 
correlations  originally  explored  by  Brown  and  Twiss.  In 
this  type  of  technique,  only  the  transverse  coherence 
length  needs  to  be  longer  than  the  spatial  separation  be¬ 
tween  detection  regions.  This  type  of  detection  can  be 
applied  to  perform  atom  antibunching  experiments  and  is 
likely  to  be  of  use  in  studying  atomic  spatial  coherence. 

The  simplicity  of  the  proposed  optical  measurement 
methods  suggests  that  experiments  based  on  these  tech¬ 
niques  may  be  useful  in  studying  quantum-mechanical 
features  of  atomic-position  measurement  in  general."  In¬ 
teresting  applications  include  the  creation  and  study  of 
one-dimensional  wave  packets,  which,  as  shown  above, 
can  be  minimum-uncertainty  Gaussian  states  under  ap¬ 
propriate  conditions.  The  evolution  of  position-squeezed 
states"  also  may  be  observable  with  one-dimensional  har¬ 
monic  wells  based  on  spatially  varying  light-shift  poten¬ 
tials. 
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APPENDIX  A.  DERIVATION  OF  THE  GREEN’S  FUNCHON  FOR  PROPAGA'HON  IN  A  UNIFORM  FORCE  FIELD 
According  to  Eq.  (14)  of  the  text,  the  required  Green’s  function  satisfies 

GF(y,yj;x,xi)=6(y  — yj)8(x  — x'j)  ,  (At) 


a  .  a  (■«  a^  .,  ,f  , 

^  ay  *  ax  2Af  dx^  ^ 


where  a  positive  infinitesimal  (e)  is  added  to  the  equation  to  ensure  causality.  The  Green’s  function  can  be  found  by 
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Fourier-transform  techniques.  We  begin  by  assuming  the  form 
G{y.y’-,x,x')=  f  f  . 

Substituting  Eq.  (A2)  into  £q.  (At)  yields 
ffdkydk^-^  ikyVy+ik^u^  +  ~k^  +  ijXo  +  € 


-if  fdkydk^ 


1  8  nk^y  +  k^x 

(2ir)^  8fcx  ^ 


'  G(k,,k^)=f  fdk^dk^-^e 


1  ikiy~y‘)^ikjx-x') 
- 7  ^ 


The  integral  involving  the  derivative  of  the  exponential,  which  arises  from  the  x-dependent  term  in  Eq.  (Al),  can  be  in¬ 
tegrated  by  parts  assuming  =  ±oo  )=0.  The  surface  terms  then  vanish,  and  the  Fourier  transform  is  found  to 

satisfy 


f  dG{ky,k^ ) 

J  aitl 


-fi  kj,v^  +  kj^v^  +  -^k^+-~XQ  —  ie  G(ky,k^)=e  '*V+*x*> 


Note  that  the  same  equation  is  obtained  by  Fourier  transforming  Eq.  (Al)  with  J  jdydx  exp  —  i{kyy+kj^x ),  provided 
thatGlj'— ►±oo,x— »  +  «)— *0.  This  equation  is  readily  solved  using  an  integrating  factor.  Let  (for  F¥=0) 

*  f  fik'^  p 

G{ky,kj,)=  A(ky,k^)exp  —‘~pf  —  >  ^^5) 


where  a  is  an  undetermined  constant.  Then, 


p  dA  {ky,k^ ) 
^  Ur 


This  has  the  solution 


ijf^  k,Oy+k:v,  +  —  +  jXo-ie  dk,  e 


fik’x  ,  F 


,  -ikyZ-ik^X- 


9t  /•''r  ft  /•^x  r 

A(ky,ky)-y  ^ ^  dk];np  i—f^  +  +  « 


-<V' 


G{ky,k„)=j,f^"dk’,'e  '*''‘exp  -ij: f ^f^y>^y+K'>x  + ~^  +  >  <A8) 

where  the  constant  a  in  Eq.  (A5)  drops  out.  Ask^  —  oo,  GIA,,*,  )-»0  for  positive  F,  since  exp(  -(ffl/F ) To 

ensure  that  G (Ac„  — ►  —  oo  )— *0,  we  must  take  6  =  —  oo .  However,  if  f  < 0,  we  need  h  =  -f  co ,  so  that  k"  >  k, .  Then, 
again  as  *±00,  G(A:^,k, )— *0. 

With  Eqs.  (A8)  and  (A2),  the  Green's  function  then  takes  the  form 

.  ..  f  r  1  ikr^y  -y'\  +  ik  X  fi  p  *x  „  ~ik"x'  .ft  f  *i  ,  ,  ,  ,  ,  ^^x  ,  F  ,,  , 

Gly,y’;x,x’)=  j  f  dkydk^—^e  dk^e  exp  +  +  +  dk,  . 


fik'y^  p 


To  incorporate  the  boundary  conditions  in  an  F-independcnt  manner,  it  is  convenient  to  change  variables,  letting 
ky={F /fi)T,  k'^=(F /fiW,  and  ky={F/fi)T".  Further,  let  b"  =  b  /[F /ft).  Note  that  if  f  >0,  dky  —  \F /ft\dT  and  if 
f  <0,  dk,^  =  —  \F/fi\dT.  In  either  case,  f  Z„dky-*\F/fi\ since  the  limits  of  integration  change  sign  with  the 
sign  of  Faccording  to  the  above  substitution.  Then,  with  P~F /ft,  the  Green’s  function  takes  the  form 


G[y,y'\x,x')=\p\  / '^’■"^'exp  -•  f^.dr'  PT'0y  + —r'^+Px^  f  df^y-^e 


1_  ikyly-y’-Vrlr-T")] 


The  last  factor  is  just  8(y  —y'  —  v,(T—T")).  To  ensure  forward  propagation,  we  want  G— *0  if  y  <y'.  This  requires 
T>  t".  Hence,  we  take  ft"  =  —  00  independent  of  the  sign  of  F.  Then  the  Green’s  function  becomes 
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G{y,y';x,x')=\p\  /”  txp-i  f'dr'  prv^  +  ^r^+Pxo  b{y -y'-v/r-r"))  . 

~  CO  ^TT  ~  CO 


By  construction,  r>  r",  so  that  the  6  function  in  Eq.  (A1 1)  is  zero  unless y  >y'.  In  this  case,  the  t"  integration  is  readi¬ 
ly  carried  out  to  obtain 

G(y,y';ac,jc')  =  e(y -y')— /”  drexp  i^TX-iplr-^ — ^  x'  exp  dr'  Pr'v^  +  ^^~ +Pxo 

u„  ~  00  I  v„  ^  y  *' ZAz 


The  r'  integration  in  the  phase  is  then  carried  out  and  yields 


G(y,y':x,x')=e(y-y')-l^  f"  dr^exp  iPiix -x')+iP^—^x'  e'*'"*  , 

-  ®  27r  Vy 


where  the  phase  4>  is  given  by 


Uj,  2  Vy 


Vy  2M  Vy  IM  Vy 


6M  v„ 


It  is  convenient  to  make  the  substitution  k  =/3r  in  Eq.  (A  1 3).  In  this  case,  for  l3 >  0,  fZ  .^dr—*  f  Z„dk /\P\.  For  P<0, 
/"^dr— — dfc/lBl  =  /!!„dfc/lj3l.  The  Green’s  function  is  then  given  by 


,y.;;,,^.,=  e(Z^/;^dk^exp  \ik  \x 


2M  Vy  2M  Vy 


Xexp  iP^-^(x'-Xo)+i-::r 

Vy  2  Vy 


6M 


The  form  of  Eq.  (A  15)  can  be  improved  somewhat  by  noting  that  in  propagating  the  final-state  amplitude,  the  integra¬ 
tion  over  x'  projects  out  the  plane-wave  component  of  wave  vector  k’,  where  k'  is  the  coefficient  of  —x'  in  the  phase 
which  appears  in  the  integrand  above.  With  the  substitution  k’  =  k  —p[{y  —y’)/Vy],  k  may  be  eliminated  in  favor  of 
the  integration  variable  k'.  Then,  eliminating  the  (')  and  using  HP^F,  one  obtains 


r-  „  1  F  l3'-:>'2ni  .iik^y-y'i 

G,{y,y,-,x,x,)  = - - - j  dk--txp  ,k  x-x,-v— - —  exp 


Vy  2M  (  Vy 


2M  u„ 


.Fy-y\.  ,  3'-3'2  1^  .  3'-y2  1' 

Xexp  (— - (Jc-Xo)  — I— -  -  ~'7T77  - 

fi  Vy  2fi  Vy  bftM  Vy 


which  is  identical  to  the  Green’s  function  given  in  Eq.  (16)  of  the  text. 


APPENDIX  B;  DERIVATION  OF  THE  INTEGRATED  RNAL-STATE  FLUX 

To  proceed,  it  is  convenient  to  perform  the  X2  integration  in  Eq.  (15)  first.  From  Eqs.  (13)  and  (15)  this  requires  eval¬ 
uation  of  the  integral 
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The  k  I  integration  required  by  Eqs.  (15),  (B3),  and  (13)  is  done  next.  This  requires  evaluation  of  the  integral 
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where  the  observation  point  y  is  chosen  downstream  from  the  Raman  region  such  thatj^  —y^  »d,  and  d  is  the  radius  of 
the  Raman  region.  Since  the  Gaussian  shape  of  the  Raman  regions  constrains  >>2  —yt  to  within  a  distance  d,  the  0 
function  which  appears  in  Gp  can  be  taken  equal  to  unity  and  theyj  integration  region  formally  extended  to  infinity. 
With  this  result  the  final-state  amplitude  takes  the  form 
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where  the  Rabi  frequencies  are  evaluated  at  the  peak  of  the  Gaussian  field  distributions. 

Our  goal  is  to  calculate  the  spatially  integrated  final-state  flux  at  large  distance  downstream  from  the  Raman  region. 
The  final-state  atomic  flux  is  proportional  to  /  =u^|  A y(y,x)P  for  the  case  of  a  collimated  atomic  beam.  The  total  sig¬ 
nal  is  then  N  ^  J"  f  dzdx  J,  where  the  integrals  extend  across  the  atomic  beam  from  —  00  to  oc.  In  order  to  take  the 
limit  y  — y,  — *  00  when  the  observation  point  is  well  outside  the  Raman  region,  it  is  necessary  to  do  the  infinite  x-z  in¬ 
tegration  before  taking  the  large-y  limit.  In  this  way,  all  of  the  flux  is  included  at  each  value  of  y.  Since  the  amplitude 
is  z  independent,  only  the  x  integral  need  be  explicitly  done.  Using  dummy  integration  variables  x\,y\  and  Jcj'.yj'  for 
the  final-state  amplitude  given  in  Eq.  (B6)  and  its  complex  conjugate,  the  required  in.egral  is  given  by 
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Note  that  the  two  terms  proportional  to  {FVy/2ii)[{y  — cancel  in  Eq.  (B6),  and  hence  do  not  appear  in  Eq. 
(B7). 

At  this  point,  the  x  integration  is  straightforward.  We  begin  by  rewriting  Eq.  (B7)  in  the  form 
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Willi  Eqs.  (B8)  and  (BOI,  one  obtains 
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In  order  to  display  explicitly  the  reality  of  the  right-hand  side  of  Eq.  (BID,  it  is  Convenient  to  make  the  substitution 
x'l  — >-x',  +(6/2).  Equation  (BID  then  takes  the  form 
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Note  that  the  right-hand  side  is  invariant  under  complex 
conjugation  if  the  substitution  y'l'^y'i  is  made  (e  changes 
sign). 

Next,  we  take  advantage  of  the  fact  that  the  observa¬ 
tion  point  y  is  chosen  downstream  from  the  Raman  re¬ 
gion  at  y,,  such  that  y  —y;  »j'i  — y,  ~^d.  In  this  case,  it 
is  convenient  to  make  the  substitution  u ,  =y2  —y,  and 
W2=j'2  — y,.  Then,  (y  — y;)"  — (y  — yj')",  which  appears 
in  the  phase  <f>  can  be  written  in  the  form 
(y  — y,  —  U| )"— (y  — yi  —  Uj)"  and  expanded  exactly  in 
powers  of  (y  — y  ,)»u,  ^d.  Hence. 
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These  results  are  substituted  into  Eq.  (B13)  and  terms  are 
collected  which  are  proportional  to  (y— y|)"  where  n 
runs  from  +2  to  —I.  It  is  found  that  all  terms  vanish 
except  for  the  one  with  n  =0.  The  phase  then  takes  the 
form 
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The  phase  does  not  depend  on  the  distance  from  the  Ra-  diffraction  and  acceleration  are  present.  Note  that  the 
man  region  y  \  to  the  observation  point  y.  This  is  reason-  phase  does  not  separate  into  two  terms  which  are  Tune¬ 
able,  since  the  integrated  flux  at  observation  points  well  tions  of  u,  and  Uj  due  to  the  cubic  dependence  on 

outside  the  Raman  region  (y— yi»d)  should  be  in-  m,— Uj.  With  Eq.  (B 15),  the  integrated  final-state  flux  is 

dependent  of  the  observation  position  even  when  given  by 
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where  e(u ,  ,U2 )  is  given  in  Eq.  (B14). 

Since  the  phase  Eq.  (BI6)  is  dependent  principally  on  the  difference  between  u,  and  it  is  convenient  to  make  the 
variable  change  u+  =(U|  +Ui)/1  and  u_  =U|  ~U2-  ^^is  case,  the  Jacobian  of  the  transformation  is  unity  so  that 
dUfdu2=du  +  du  Integration  is  over  the  entire  plane.  Substituting  u ,  =u  + -t-(u  _ /2)  and  U2  =  u+-(u_/2)  into 
Eq.  (B16),  one  obtains 
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Equation  (BI7)  gives  the  general  result  for  the  integrated  flux  in  the  final  state  in  terms  of  the  initial-state  amplitude. 
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